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University  of  Washington 
Abstract 

Use  of  Power  Balance  to  Model  Converging  Flows  in  Bilateral  Junctions  for  Ventilation 

Systems 

by  James  P.  Curran 

Chairperson  of  the  Supervisory  Committee:  Professor  Steven  E.  Guffey 

Department  of  Environmental  Health 

Little  information  is  available  on  the  nature  of  flow  in  bilateral  junctions  for  industrial 
ventilation  systems.  The  American  Conference  of  Governmental  Industrial  Hygienists 
recommends  against  bilateral  junctions  without  providing  any  justification.  The  American 
Society  of  Heating,  Refrigeration  and  Air-conditioning  Engineers  provides  a  table  of 
coefficients  for  one  type  of  bilateral  junction;  this  table  was  found  to  be  in  error. 

The  main  objective  of  this  study  was  to  determine  if  the  power  balance  method  and  the 
linear  kinetic  power  model  could  be  used  to  characterize  and  predict  the  flows  in  bilateral 
junctions.  A  secondary  goal  was  to  determine  the  relative  energy  efficiency  of  bilateral 
junctions  compared  to  two  single  lateral  junctions  in  series. 

Power  losses  in  one  bilateral  junction  were  experimentally  determined  for  35  conditions 
where  the  velocity  ratios  of  the  converging  ducts  were  varied  from  0-3.7.  The  linear  kinetic 
power  model  provided  a  good  fit  to  the  empirical  data  (Adjusted  =  90.9%).  An  adjusted 
r2  of  97.7%  was  obtained  for  an  important  subset  of  data  (velocity  ratios  ranging  from 
0.75  to  1.25).  As  expected,  power  losses  were  lowest  when  velocity  ratios  were  near  one. 
The  power  losses  in  the  bilateral  junction  tested  were  45-92%  higher  than  the  predicted 
power  losses  for  two  single  lateral  junctions  in  series.  This  increase  in  power  loss  will  have 
an  insignificant  impact  on  the  overall  pressure  requirements  for  most  ventilation  systems. 

In  this  study,  the  maximum  increase  in  total  pressure  due  to  the  bilateral  junction  was 
estimated  as  3.5%  higher  than  the  predicted  total  pressure  in  a  hypothetical  setup  containing 
two  single-lateral  junctions. 
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CHAPTER  I 
INTRODUCTION 


Background 

Currently  most  industrial  exhaust  ventilation  systems  are  designed  using  one  of  two 
methods  to  estimate  system  pressure  requirements.  In  the  American  Conference  of 
Governmental  Industrial  Hygienist's  (ACGIH)  method,  "static  pressure  losses"  are 
calculated  for  system  components.^)  "Total  pressure  losses”  are  calculated  in  the  American 
Society  of  Heating  and  Air-conditioning  Engineer's  (ASHRAE)  method.^^)  These  methods 
have  much  in  common  since  total  pressure  (TP)  can  be  calculated  by  adding  the  velocity 
pressure  (VP)  to  the  static  pressure  (SP)  and  since  they  both  compute  losses  in  terms  of 
velocity  pressure  coefficients.  On  the  other  hand,  these  methods  differ  in  the  way  they 
estimate  the  distribution  of  flows  and  assign  "pressure  losses"  at  junctions.  The  ACGIH 
method  assumes  that  flows  are  determined  by  a  "governing"  static  pressure  in  the  junction 
(SPj).  while  the  ASHRAE  method  assumes  the  total  pressure  downstream  from  a  junction 
(TPm)  must  be  the  same  for  all  paths  through  the  junction.  Both  methods  assume  junction 
pressure  losses  are  proportional  to  velocity  pressures. 

The  power  balance  method  and  the  linear  kinetic  power  model  provide  a  third  method 
of  estimating  the  distribution  of  flows  and  the  losses  in  a  junction.^^-^)  The  power  balance 
model  is  an  application  of  the  law  of  conservation  of  energy  to  determine  the  energy  change 
per  unit  time  (i.e.  power)  for  conditions  typical  in  converging  flow  junctions. (3)  The  linear 
kinetic  power  model  assumes  that  the  power  losses  in  any  junction  are  related  linearly  with 
the  kinetic  powers  of  each  branch  entering  the  junction. 

ACGIH  currently  does  not  provide  information  to  predict  flows  and  pressures  in 
bilateral  junctions  (see  Figure  1)  and  recomn:ends  against  them  without  providing  any 
justification.  ASHRAE  provides  information  on  a  very  limited  subset  bilateral  junction 
configurations.  The  kinetic  power  model  should  theoretically  apply  to  bilateral  junctions  but 
only  single-lateral  junctions  (see  Figure  2)  have  been  tested  with  this  method. 

There  are  many  practical  applications  for  bilateral  junctions;  but,  in  the  absence  of  data 
for  design  their  use  is  problematic.  Bilateral  junctions  may  provide  economic  benefits; 


however,  the  potential  costs  of  using  these  junctions  (i.e.  increased  pressure  requirements, 
increased  settling,  etc.)  are  not  currently  known. 


Research  Objectives 

The  primary  goal  of  this  research  was  to  characterize  the  flow  in  a  bilateral  Junction  and 
to  determine  if  the  kinetic  power  model  could  be  successfully  applied  to  predict  power 
losses.  If  successful,  similar  research  could  be  completed  on  many  bilateral  junctions  with 
various  geometric  variables  to  develop  a  general  model  applying  to  a  wide  range  of  bilateral 
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(4) 

-  junctions  as  has  been  done  with  single-lateral  junctions.  ’  A  secondary  goal  was  to  test  the 
hypothesis  that  the  power  losses  for  three  flows  joined  at  a  bilateral  junction  are  not 
significantly  different  than  the  combined  power  losses  of  identical  flows  joined  by  two 
single-lateral  junctions. 


Current  Design  Methods 

ACGTH  Method 

In  Industrial  Ventilation,  the  ACGIH  provides  a  method  for  designing  systems 
containing  single  lateral  junctions.(^)  ACGIH  recommends  against  bilateral  junctions,  but 
provides  no  justification  for  this  recommendation  and  does  not  provide  a  method  of 
designing  systems  containing  these  junctions.^^^  However,  ACGIH  does  provide  two 
alternatives  to  bilateral  junctions.  The  "preferred  option"  is  to  place  two  single  lateral 
junctions  in  series  very  close  together  as  shown  in  Figure  3.  The  second  "acceptable 
option"  (Figure  4)  is  similar  to  a  bilateral  junction  except  that  the  side  laterals  are  offset 
from  each  other.  Appendix  A  contains  the  original  figure  from  the  ACGIH  manual 
showing  these  recommendations. 


ASHRAE  Method 

ASHRAE  provides  loss  coefficients  for  a  single  type  of  bilateral  junction  in  its  1989 
Fundamentals  Handbook:  Appendix  B  contains  a  copy  of  the  table  of  coefficients. (2)  In 
the  ASHRAE  bilateral  junctions,  the  diameters  of  laterals  are  equal  to  each  other  and  the 
coaxial  ducts  also  have  equal  diameters.  This  type  of  junction  may  be  useful  for  general 
room  ventilation  in  a  plenum  type  system;  but  its  value  is  limited  in  industrial  ventilation 
design  particularly  in  systems  carrying  particulates.  In  these  systems,  maintaining 
minimum  target  velocities  is  critical  to  prevent  settling.  The  bilateral  junction  cited  by 
ASHRAE  forces  designers  toward  systems  with  excessive  velocities  in  the  duct 
downstream  of  the  junction  since  the  combined  areas  of  the  upstream  ducts  are  greater 
than  the  downstream  duct.  These  excessive  velocities  increase  system  operating  and 


Junctions 


Junctions 


maintenance  costs  because  of  higher  power  requirements  and  abrasion.  In  gas  and  vapor 
conveying  systems,  minimum  target  velocities  are  not  critical;  but  similar  tradeoffs  occur 
between  initial  costs  and  long-term  operating  costs. 
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ASHRAE  adapted  their  coefficients  from  Idelchik.^^^  When  transferring  Idelchik's  table 
c-t  vcfficients,  ASHRAE  made  an  error  in  the  coefficients  for  the  "main"  portion  of  the 
table  by  printing  the  ratio  of  the  areas  for  the  side  branches  (A2b/Aj5)  as  an  independent 
variable.  Idelchik's  original  table  lists  the  ratio  of  flows  (Q2b/Qlb)  the  independent 
variable.  Idelchik  provides  several  tables  of  coefficients  for  similar  junctions  with  different 
branch  entry  angles  (see  Appendix  C). 

The  author  verified  Idelchik's  table  of  coefficients  that  was  adapted  by  ASHRAE  by 
substituting  the  independent  variables  from  the  tables  into  the  equations  accompanying  the 
tables.  Several  differences  were  found  between  the  calculated  loss  coefficients  and  the  table 
values.  In  most  cases,  the  difference  between  the  table  and  calculated  values  was  only 
0.01-0.03  (dimensionless  coefficient);  however,  in  several  cases  the  difference  was 
between  0. 10  and  0.50.  In  one  case  the  table  value  was  10  times  the  calculated  value. 
Appendix  D  contains  the  tables  of  corrected  coefficients.  No  effort  was  made  to  verify  the 
coefficients  for  any  of  the  other  bilateral  junctions  provided  by  Idelchik. 


Power  Balance  Method 

Air  flowing  through  a  ventilation  system  has  three  energy  components:  kinetic, 
potential  and  intemal.(^)  The  static  pressure  between  two  points  in  a  system  provides  the 
driving  force  for  air  and  is  thus  the  source  of  potential  energy.  The  air  moves  from  one 
point  to  another  at  a  certain  velocity  and  thus  has  kinetic  energy.  Finally,  air  has  internal 
energy  that  changes  because  of  heat  transfer.  At  any  cross-section,  the  air  flows  through  a 
system  at  a  certain  rate.  As  a  result,  these  energies  can  be  expressed  as  an  energy  rate,  or 
power. 

The  general  form  of  the  total  power  of  flow  at  a  cross- section  (Equation  1)  is  quite 
complicated;  however,  assumptions  can  be  made  which  greatly  simplify  the  general 
equation  without  significant  loss  in  accuracy  for  typical  ventilation  systems.(^)  Following 
is  the  general  equation  for  the  total  power  flow  at  a  cross  section:(^ 
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TFP=  J.(vxp)  [  — +  e  +  5(vx^  +  Vy^  +  V22)  +  hzg]dA  . (1) 

A  p 

where:  TFP  =  total  flow  power  at  a  cross  section  of  duct 

Vx  =  velocity  of  air  along  the  axis  of  the  duct 
p  =  mass  density  of  air 

S  P  =  static  pressure  at  the  cross  section 

e  =  specific  internal  energy  of  the  airstream 

Vy  =  air  velocity  toward  the  side  of  the  duct  along  the  y 
axis 

v^  =  air  velocity  toward  the  side  of  the  duct  along  the  z 
axis 

=  elevation  of  the  flow  from  a  datum 
g  =  acceleration  due  to  gravity 

For  any  junction,  where  no  external  power  is  added,  the  power  of  the  airflow  exiting 
the  junction  must  equal  the  sum  of  the  powers  of  the  incoming  airflows  minus  the  power 
lost  (dissipated  as  heat).  Power  may  also  be  lost  from  the  system  by  the  conduction  of  heat 
through  the  walls  of  the  junction  to  the  surrounding  environment.  Equation  1  may  be 
simplified  based  on  the  following  assumptions:(^) 

1.  One  dimensional  flow.  Air  flows  axially  along  the  duct;  therefore,  Vy  and  v^ 
are  zero. 

2.  Incompressible  flow.  Negligible  errors  result  from  this  assumption  for  velocities 
up  to  Mach  0.2  Mach  0.2  is  equivalent  to  13,200  feet  per  minute.  The 
maximum  velocity  used  in  this  study  was  less  than  one  half  this  value,  yet 
represents  an  extreme  value  for  practical  ventilations. 

3.  Adiabatic  conditions.  Since  the  air  in  the  duct  comes  from  within  the  laboratory, 
there  is  very  little  difference  (±2°C)  between  the  temperature  of  the  air  in  the  duct 
and  the  environment  surrounding  the  duct.  As  a  result,  there  is  no  large 
temperatiu’e  gradient  to  effect  a  significant  amount  of  conduction  across  the 
junction  or  across  the  ducts  between  the  junction  and  the  measuring  stations.  In 
addition,  the  lengths  of  duct  are  relatively  short  which  minimized  the  opportunity 
for  heat  transfer. 


4.  Uniform  flow.  This  assumption  can  be  verified  by  observing  the  uniformity  of  the 
velocity  profiles  during  data  collection.  If  the  velocity  profiles  are  reasonably 
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uniform,  then  the  air  density  and  static  pressure  should  also  be  uniform.  In  fact,  it 
is  not  necessary  for  the  velocity  profiles  to  be  uniform  as  long  as  the  velocity 
contours  are  nearly  the  same  upstream  and  downstream.(3) 

5.  No  change  in  elevation.  Since  the  density  of  air  is  low,  negligible  errors  will 
result  from  small  elevation  changes  within  the  laboratory.  All  measurement 
stations  were  within  1  inch  of  the  same  elevation  in  this  experiment. 


Qb 


Figure  5:  Bilateral  Junction  with  Nomenclature 


Following  is  the  simplified  Power  Balance  equation  for  a  bilateral  junction  (see  Figure 

5):® 

LP  =  Qa(SPa+VPa)  +  ObCSPb+VP^)  +  OsiSPs+VP^)  -  Qm(SPm+VPm)  . (2) 

where:  LP  =  total  lost  power 
Q  =  flow  in  the  duct 
SP  =  static  pressure  in  the  duct 
VP  =  average  velocity  pressure  in  the  duct 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in  ducts  A,  B, 
S  and  M 

Rows  and  pressures  cannot  be  measured  near  the  junction  because  of  the  disturbances 
induced  by  the  colliding  airstreams.  As  a  result,  these  measurements  must  be  taken  at  some 
distance  from  the  junction.  Measurements  taken  far  upstream  and  downstream  include  the 


losses  attributable  not  only  to  the  junction  but  to  friction  as  well.  To  find  the  power  loss 
attributable  to  the  Junction  alone,  the  power  lost  due  to  friction  between  the  points  where 
pressure  measurements  are  taken  and  the  junction  must  be  estimated  and  then  subtracted 
from  the  total  dissipated  power.  The  expression  for  the  power  lost  due  to  friction  follows: 

LPF  =  QmfASPFj.n,)  +  QaCASPFj.^)  +  Qb(ASPFj.b)  +  QsfASPFj.^)  . (3) 

where:  LPF  =  power  lost  due  to  friction 

ASPFj.^  =  predicted  change  in  static  pressure  (due  to  friction) 
between  the  junction  and  the  duct  measuring  station 

Q  =  airflow  in  the  duct 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in  ducts 
A,B,  S  and  M 

Finally,  the  lost  power  attributable  to  a  junction  can  be  determined  empirically  from  the 
following: 


LPJ  =  LP-LPF  . (4) 

where:  LPJ  =  power  lost  due  to  the  junction 
LPF  =  power  lost  due  to  friction 
LP  =  total  lost  power 


By  expansion,  Equation  4  becomes: 

LPJ  =  [Qa(SPa+VPa)  -h  Qb(SPb+VPb)  +  Qs(SPs+VP,)  -Q^CSP^+VP^)]  . (5) 

-  [QmfASPFj.n,)  Qa(ASPFj.a)  -h  Qb(ASPFj.b)  +  QsCASPFj.^)] 

where:  LPJ  =  power  lost  due  to  the  junction 

ASPFj.j  =  predicted  change  in  static  pressure  (due  to  friction) 
between  the  junction  and  the  duct  measuring  station 


Q  =  airflow  in  the  duct 
SP  =  static  pressure  in  the  duct 
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VP  =  average  velocity  pressure  in  the  duct 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in  ducts 
A,  B,  S  and  M 

Guffey  and  Fraser  determined  that  the  power  losses  in  a  single  lateral  junction  were  a 
linear  function  of  the  kinetic  powers  of  the  incoming  and  outgoing  ducts.(^)  Applying  the 
kinetic  power  model  to  a  bilateral  junction,  the  following  model  results; 

LPJ  =  QaVPaKa  +  QbVPfaKb  +  QsVPgKs  . (6) 

where:  LPJ  =  power  dissipated  across  the  junction 

Q  =  air  flow  through  the  duct 

VP  =  average  velocity  pressure  in  the  duct 

K  is  an  empirical  loss  coefficient  for  the  duct 

Subscripts  a,  b  and  represent  reference  cross-sections  in  ducts  A,  B 
and  S 

The  Case  for  Bilateral  Junctions 

Bilateral  junction  systems  have  practical  advantages  over  multiple  single-lateral 
systems.  Normally,  plant  equipment  is  located  to  provide  the  most  efficient  flow  for  the 
work  process  and  to  conserve  floor  space.  As  a  result,  equipment  is  often  positioned  on  the 
plant  floor  symmetrically.  If  plant  processes  require  local  exhaust  ventilation  in  this 
situation,  a  bilateral  junction  to  merge  system  branches  is  an  attractive  alternative  since  it 
simplifies  system  design  and  installation  (see  Figure  6). 


Figure  6:  Bilateral  Junction  Design  Option 
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In  reality,  the  ACGIH  alternatives  to  bilateral  junctions  (Figures  3  and  4)  may  be  no 
more  efficient  than  bilateral  junctions.  One  might  assume  that  the  loss  coefficients  for  single 
lateral  junctions  apply  to  these  alternatives  since  ACGIH  provides  no  correction  factors  or 
further  explanations.  However,  this  is  a  precarious  assumption  since  locating  system 
components  too  close  together  can  cause  calculated  pressure  losses  to  be  in  error  by  a  factor 
of  two  or  more  because  of  interaction  effects.(^)  Miller  suggests  an  interaction  correction 
factor  of  unity  can  be  used  only  when  the  junctions  are  three  or  more  diameters  apart.(*) 

The  ACGIH  "preferred"  option  does  not  appear  to  meet  this  guideline  (see  Appendix  A). 

Since  the  behavior  of  flow  in  bilateral  junctions  has  not  been  well  characterized,  and 
since  they  are  not  recommended  by  the  ACGIH,  designers  are  forced  into  less  economical 
alternate  designs.  Worse  yet,  designers  may  be  forced  into  guessing  what  will  happen  if 
these  junctions  are  used  and  then  to  adjusting  the  system  after  construction  to  compensate 
for  pressure  loss  calculation  errors.  Figures  7  and  8  show  two  alternate  system  designs  that 
avoid  bilateral  junctions.  The  parallel  systems  design  in  Figure  7  has  the  disadvantages  of 
higher  initial  costs  and  additional  space  requirements  that  may  lead  to  interferences  with 
other  plant  equipment  or  utilities.  The  offset  single-lateral  junction  design  shown  in  Figure 
8  is  more  economical  and  consumes  less  space  than  the  parallel  system  design;  however,  it 
IS  still  inferior  to  the  bilateral  junction  design  with  regards  to  initial  costs  and  space 
requirements.  Long-term  operating  costs  will  also  be  slightly  higher  in  the  offset  single¬ 
lateral  system  than  the  bilateral  system  if  the  longer  laterals  increase  system  pressure 
requirements.  A  final  factor  that  must  be  considered  with  the  offset  single-lateral  design  is 
the  unknown  increase  in  pressure  losses  that  may  result  because  of  interaction  effects  if  the 
single-lateral  junctions  are  located  close  together. 


Figure?:  Parallel  Systems  Design  Option 


CHAPTER  II 

MATERIALS  AND  METHODS 


Laboratory  Setup 

The  apparatus  included  the  fan;  the  ducts  and  their  connectors;  the  supporting  frame  for 
the  ducts;  the  velocity  pressure  measurement  stations;  Pitot  tubes  and  associated  tubing; 
pressure  sensors;  and  the  computer  and  its  peripheral  devices  for  data  collection. 


Fan,  Motor  and  Transmission 

Airflow  was  provided  by  a  centrifugal  fan  (Buffalo  Forge,  MW40,  Buffalo  NY)  and 
was  varied  using  an  infinitely  variable  geared  transmission  device  (Speed  Control 
Industries  Inc,  Model  7-lOA  (E)  0:1-1  CCW,  Richland  WA)  connected  to  a  20  horsepower 
motor.  The  drive  shaft  of  the  fan  was  connected  to  the  transmission  by  belts  and  pulleys, 
and  the  shafts  of  the  motor  and  transmission  were  connected  shaft  to  shaft  with  a  locking 
device.  Fan  rotation  rates  were  varied  from  about  200  to  16(K)  revolutions  per  minute 
(rpm).  Rotation  rates  drifted  from  a  nominally  fixed  rate  by  an  average  of  7.5  rpm/hour 
(range  2-28  rpm/hr  over  a  30  minute  period).  In  general,  the  run-time  per  observation 
period  was  approximately  30  minutes  for  friction  and  2  hours  for  junctions.  The  drift 
tended  to  be  highest  at  higher  rotation  rates. 


Duct  System 

The  ducts  used  in  this  study  were  standard  20  and  24  gauge  spiral  wound  galvanized 
steel  duct  that  was  purchased  from  a  manufacturer  of  ventilation  system  components 
(United  McGill,  Stockton  CA).  All  ducts  were  cleaned  with  hot  water  and  dishwashing 
detergent  to  remove  a  vegetable  oil  coating  applied  during  the  manufacturing  process. 

Ten  foot  lengths  (±  3  inches)  of  spiral  wound  duct  rested  on  0.75  inch  cabinet-grade 
plywood  vertical  spacers  (Figure  9)  that  were  mounted  every  10  feet  on  an  80  foot  long 
metal  frame  suspended  30  inches  from  the  laboratory  ceiling.  The  frame  was  supported 
every  6  feet  with  0.375  inch  steel  threaded  rod,  and  was  leveled  with  a  rotating  laser.  Each 
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vertical  spacer  was  centered  on  the  longitudinal  axis  of  the  support  frame.  Separate  hangers 
spaced  at  6-8  foot  intervals  were  installed  to  support  the  system  branches  (see  Figure  10). 


Figure  9:  Photograph  of  Plywood  Vertical  Spacers 


The  ducts  were  connected  to  each  other  with  18  gauge  galvanized  steel  insert  couplers 
manufactured  for  that  purpose.  The  couplers  were  4  inches  long,  and  v/cre  approximately 
equal  to  the  nominal  duct  diameter  at  the  midpoint  and  0. 1 8  inches  less  than  the  nominal 
duct  diameter  at  each  end  to  allow  a  tight-fitting  insertion.  The  ends  of  the  ducts  connected 
to  the  measurement  stations  terminated  in  2  inch  iron  flanges  constructed  from  0.14  inch 
thick  angle  iron.  The  upstream  end  of  the  each  velocity  pressure  measurement  station 
terminated  with  a  plywood  (0.75  inch  thick)  over-sized  rectangular  flange.  The  ducts  were 
aligned  to  the  measurement  station  by  bolting  each  duct  flange  to  a  corresponding 
measurement  station  flange.  No  part  of  the  flanged  connectors  projected  into  the  airflow. 
These  connectors  provided  very  good  alignment  and  minimized  potential  disturbances  close 
to  the  measuring  station.  Figure  1 1  shows  the  standard  connector  on  the  left  and  the 
flanged  connector  on  the  right. 


All  duct-to-duct  connections  were  sealed  with  duct  tape  to  prevent  leaks.  The  duct-to- 
measurement  station  connections  were  sealed  with  a  combination  of  duct  tape  and  silicone 
sealant.  Ducts  were  aligned  in  the  horizontal  plane  using  a  string  line  and  their  alignment  in 
the  vertical  plane  was  verified  with  a  construction-type  level. 

Pressure  Measurement  Stations 

A  "drawn-over-mandrel"  (DOM)  0.25  or  0.5  inch  thick  steel  tube  with  an  inner 
diameter  matching  the  nominal  inner  diameters  of  the  duct  was  used  for  all  velocity 
pressure  traverses.  All  static  pressures  and  static  pressure  differentials  were  measured  at  the 
duct  centerline. 

All  pressure  measurements  were  taken  using  ANSI/ASHRAE  Standard  41.2-1987 
hemispherical  head  0. 125  inch  Pitot  rubes  mounted  in  Pitot  tube  holders  with  log-linear 
insertion  depths  controlled  to  within  0.0005  inches  (see  Figures  12  and  13). (9, 10,1 1,12)  Xq 
minimize  leaks  around  the  Pitot  tubes,  a  small  (0.5  inch  by  1.0  inch)  galvanized  steel  sheet 
metal  strip  with  a  0. 15  inch  diameter  hole  drilled  in  it  was  taped  over  the  hole  that  had  been 
drilled  in  the  duct  or  DOM  tubing  before  inserting  the  Pitot  tube.  All  velocity  pressure 
measuring  stations  contained  three  Pitot  tubes  for  conducting  10-point  traverses  along  3 
diameters  off-set  by  120  degrees  but  in  the  same  plane  (see  Figure  14).  Each  of  the  Pitot 
tube  holders  was  mounted  to  the  DOM  tubing  as  shown  in  Figure  12. 


Leak  tests 

Pressure  was  conducted  from  the  Pitot  tubes  to  the  pressure  sensors  using  3/16  inch 
inside  diameter  (1/16  inch  wall  thickness)  Tygon®  tubing  and  low  density  polypropylene 
quick-disconnect  fittings.  The  number  of  fittings  was  kept  to  a  minimum  to  reduce  the 
potential  for  leaks.  Each  pathway  was  tested  in  place  for  leaks  with  approximately  4  inches 
water  gage  (w.g.)  pressure.  Over  a  5  minute  span,  none  showed  measurable  pressure 
losses  (<0.005  inches  w.g.). 


Back  View  Side  View  Front  View 


Figure  12:  Diagram  of  Pitot  Tube  Holding  Device  and  Mounting 


Figure  13:  Photograph  of  Pitot  Tube  Holding  Device 


17 


DOM  tubing 


3  Pitot  tubes 
side  view 


A 


end  view 


120  deg 


Figure  14:  Velocity  Pressure  Measuring  Station  Setup 


A  leak  test  of  ducts  was  completed  on  each  setup  before  collecting  data  by  sealing  all 
duct  inlets  with  duct  tape  and  developing  a  static  pressure  of  approximately  2.5  inches 
within  the  system.  Leaks  were  sought  and  sealed  until  the  differential  static  pressures 
across  the  experimental  lengths  were  less  than  0.005  inches  and  there  was  no  measurable 
centerline  velocity  pressure.  After  completing  the  leak  test,  the  duct  tape  was  removed  from 
the  end  of  the  ducts  and  a  steel  ring  with  open-celled  honeycomb  air-flow  straightener  was 
placed  on  the  end  of  each  duct  not  having  this  material  already  installed  several  diameters 
downstream  of  the  opening. 


Computer  Data  Acqui.sition  System  and  Data  Reduction 

Pressure  transducers  were  connected  to  a  computer  (Apple  Macintosh™  SE)  containing 
a  data  acquisition  card  via  a  terminal  panel  (Strawberry  Tree,  Sunnyvale  CA).  The 
acquisition  card  converted  analog  voltages  to  digital  inputs  that  were  recorded  to  a  computer 
file  using  "Analog  Connection  Workbench™  2.0"  software  (Strawberry  Tree,  Sunnyvale 
CA).  A  binary-coded  thumb  wheel  (Cherry  Electrical  Products  Corporation,  Waukegan  IL) 
was  used  to  identify  groups  of  data  and  individual  data  points.  The  Workbench™ 
computer  program  produced  ASCII  text  files  containing  the  date,  time,  identification  code 
and  transducer  voltage  output  for  each  of  eight  channels.  Figure  15  is  a  wiring  diagram  for 
the  transducers  and  Figure  16  is  a  wiring  diagram  for  the  transducers. 


7.6v  source 


During  transducer  calibration,  individual  calibration  points  were  consecutively 
numbered  with  settings  on  the  thumb  wheel  ranging  from  0  through  79.  Normally  no  more 
than  50  calibration  points  were  accomplished  in  any  one  day. 

During  friction  and  junction  data  collection,  groups  of  data  were  identified  by  setting 
the  thumbwheel  to  a  "run"  number  from  41  through  75  and  pressing  the  "pulse  on/off ' 
switch  while  the  switches  for  all  transducers  were  in  the  off  position.  The  thumbwheel  was 
then  set  to  a  number  (1-39)  to  identify  an  individual  data  point  within  a  given  run.  Only  one 
analog  channel  was  switched  on  at  a  time.  Table  1  contains  the  coding  system  used  to 
identify  data  "runs"  and  individual  data  points.  When  all  data  for  a  "run"  was  collected  the 
signal  switches  for  all  transducers  were  turned  off  and  the  next  run  number  was  logged. 


Table  1:  Computer  Acquired  Data  Coding  System 


Pressure  Type 

Velocity  pressure  traverse  points 

1-30 

Centerline  velocity  pressures 

31-33 

Static  pressure  at  VP  measuring  location 

34 

Differential  static  pressures 

35-39 

Calibration  pressures 

76-79 

Computer- generated  ASCII  text  files  were  analyzed  using  one  of  two  computer 
programs  developed  by  Dr.  Guffey.  Calibration  data  was  analyzed  using  the  "Calibration 
Digestion"  program.  Friction  and  junction  data  were  analyzed  using  the  "Digest  Workbench 
Data"  program. 

The  "Digest  Workbench  Data  Program"  reduced  data  to  the  average  velocity,  average 
velocity  pressure,  and  differential  static  pressures.  The  program  opened  and  read  the  ASCII 
text  file  containing  transducer  voltages;  converted  all  voltages  to  pressures  based  on 
calibration  data;  opened  and  read  the  text  file  containing  the  laboratory  dry  bulb  and  wet 
bulb  temperatures,  barometric  pressure,  and  airstream  temperature;  calculated  the  corrected 
barometric  pressure;  calculated  the  airstream  density  factor  from  the  duct  temperature,  water 
vapor  concentration,  duct  static  pressure  and  corrected  barometric  pressure;  and  calculated 
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the  average  velocity  and  average  velocity  pressure  for  each  traverse.  The  average  velocity 
pressure  in  this  instance  is  not  the  average  of  the  velocity  pressures,  but  the  average 
velocity  pressure  calculated  from  the  average  velocity  in  the  duct. 

While  the  "Digest  Workbench  Data"  program  was  running,  it  graphed  velocity  contours 
for  each  traverse  on  the  computer  screen,  allowing  identification  of  possible  data  collection 
errors  and  poor  flow  conditions.  Figure  17  is  a  photograph  of  the  screen  output  during  data 
analyses.  The  two  traverses  on  the  right  side  of  Figure  17  have  uniform  velocity  profiles 
while  the  left  side  shows  a  case  where  a  leak  developed  at  the  tubing/Pitot  tube  interface. 

The  program  produced  three  text  files.  One  file  contained  information  on  each 
individual  velocity  pressure;  this  file  could  be  used  to  reexamine  velocity  profiles.  The 
second  file  contained  the  following  information  for  each  run: 

Date 

Time 

Run  number 

Airstream  temperature 

Corrected  barometric  pressure 

Room  density  factor 

Number  of  traverse  points 

Duct  density  factor 

Average  velocity  pressure 

Average  centerline  velocity  pressure 

Each  centerline  velocity  pressure 

"True"  velocity  pressure 

Static  pressure  at  the  VP  measuring  station 

All  differential  static  pressures  measured 

The  third  file  was  a  self  calibration  file  that  contained  recorded  voltages  and  predicted 
pressures.  This  file  was  used  to  develop  calibration  curves  for  each  data  collection  set  using 
data  points  where  both  inclined  manometer  and  a  transducer  measurements  were  taken. 

This  procedure  is  discussed  later  in  the  equipment  calibration  portion  of  the  methods. 
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The  "Digest  Workbench  Data"  program  eliminated  calculation  errors  that  could  be 
expected  from  reducing  this  large  quantity  of  data  by  hand.  The  program  also  provided  a 
better  method  of  auditing  calculations  than  is  available  with  computer  spreadsheet 
programs. 


Figure  17:  Velocity  Contour  Output  from  Digest  Workbench  Data  Program 


Measuring  Equipment  and  Calibration 

Pressure  Measuring  Equipment 

Before  pressure  data  were  collected,  pressure  measuring  devices  were  calibrated  with  a 
hook  gage  (Dwyer  Instruments  Inc,  Series  1425,  Michigan  City  IN)  having  a  resolution  of 
0.(X)1  inches  w.g.  The  hook  gage  contained  fluorescein  green  color  concentrate  that  was 
mixed  with  distilled  water  according  to  the  manufacturer's  instructions.(^^)  Before  each 
use,  the  hook  gage  was  leveled  using  the  built  in  bubble  levels  and  the  hooks  were  zeroed. 

Two  inclined  manometers  (Meriam  Instrument,  Model  40HE35WM,  Cleveland  OH) 
with  resolutions  of  0.005  inches  w.g.  were  used  as  the  predominant  measuring  devices  for 
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static  and  differential  pressures  4  inches  w.g.  or  less.  Each  manometer  contained  Meriam 
1000  Green  Concentrate  (No.  922WA)  mixed  with  distilled  water  according  to  the 
manufacturer's  instructions.^'"^)  The  inclined  manometers  were  mounted  on  a  wall  and 
leveled  using  the  built-in  bubble  level.  The  level  and  zero  for  each  inclined  manometer  was 
checked  and  adjusted  as  necessary  before  each  set  of  readings. 

A  Meriam  Instruments  Meri-cal®  Model  DP2(X)I  produced  a  digital  readout  with  a 
display  resolution  of  0.01  inches  w.g.  It  was  electronically  zeroed  with  the  zero  switch 
before  each  use  and  was  used  with  internal  damping  toggled  on.  The  Meri-cal®  was  very 
sensitive  to  pressure  changes  and  occasionally  displayed  wide  fluctuations  in  the  digital 
readout  even  with  the  damping  switch  on.  An  attempt  to  provide  physical  damping  by 
placing  charcoal  tubes  in  the  Tygon®  tubing  leading  from  the  Meri-cal®  appeared  to  have 
little  effect.  Because  of  its  low  resolution  and  readout  fluctuations,  the  Meri-cal®  was  used 
only  to  provide  secondary  measurements  for  pressures  exceeding  4  inches  w.g.  When  the 
Meri-cal®  was  used,  the  mean  of  the  high  and  low  readouts  was  recorded  as  the  pressure. 

Two  models  of  pressure  transducers  (Omega  Engineering  Inc.,  PX160  and  PX154 
Series,  Stamford  CT)  were  used  as  the  primary  velocity  pressure  measuring  devices  in  this 
study.  The  transducers  were  mounted  on  a  plywood  panel  that  was  affixed  to  the  laboratory 
wall  (see  Figure  18).  The  transducers  produced  voltages  or  currents  proportional  to  applied 
pressures.  The  signal  output  was  digitized  and  recorded  to  a  computer  file  by  the  data 
acquisition  system.  Transducers  were  allowed  to  warm-up  at  least  one-half  hour  before 
each  use. 

The  hook  gage,  hand  pump  (Meriam  Instruments,  Model  B34348),  and  all  pressure 
sensing  instruments,  were  connected  to  a  glass  manifold  with  Tygon®  tubing  for 
calibration.  As  a  result,  each  instrument  experienced  the  same  pressure  developed  by  the 
hand  pump.  All  connections  were  sealed  with  a  silicone  sealant  and  the  calibration  setup 
was  tested  for  leaks  at  the  start  of  each  run  using  approximately  3  to  4  inches  w.g.  of 
positive  pressure.  The  setup  was  assumed  to  be  free  of  significant  leaks  if  the  leakage  rate 
was  less  than  0.0001  inches  w.g.  per  minute.  This  assumption  is  based  on  the  fact  that  all 
instruments  were  read  for  each  calibration  point  in  a  period  of  2  minutes  or  less. 
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Figure  18:  Photograph  of  the  Pressure  Transducers  Used 


The  room  temperature  and  barometric  pressure  were  recorded  at  the  start,  the  finish, 
and  periodically  during  each  calibration  run.  Generally,  calibration  runs  started  at  the  low 
end  of  the  hook  gage  range  and  progressed  to  the  high  end  of  the  range  in  one  tenth  of  an 
inch  w.g.  increments.  Once  the  high  end  of  the  range  was  reached,  additional  calibration 
points  were  obtained  by  successively  releasing  pressure  from  the  system  until  the  original 
pressure  was  achieved.  The  hook  gage  pressure  and  vacuum  readings,  inclined  manometer 
readings,  and  the  Meri-cal®  readings  were  recorded  on  a  data  collection  worksheet  for 
each  calibration  point  (see  Appendix  E).  The  hook  gage  range,  which  depended  on  the 
spacer  rod  selected,  was  also  noted  for  each  run.  The  voltage  produced  by  each  pressure 
transducer  for  each  calibration  point  was  recorded  to  a  computer  file  using  the  computer 
data  acquisition  apparatus  and  software.  To  avoid  the  potential  for  electronic  interference 
between  transducers,  the  voltages  resulting  from  the  applied  pressures  were  recorded  for 
each  transducer  with  the  signal  from  all  other  transducers  switched  off. 


Numerous  calibration  runs  ever  a  period  of  five  months  were  completed  before  data 
collection  began.  Many  calibration  runs  were  conducted  without  the  Meri-cal®  because  it 


was  out  for  repair.  Table  2  provides  a  summary  of  the  number  of  calibration  points  for 
each  piece  of  equipment  used  in  this  study. 
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The  hook  gage  pressure  values  were  corrected  for  the  slight  changes  in  indicating  fluid 
density  due  to  temperature  variation.  Figure  19  shows  the  relationship  between  the  hook 
gage  correction  factor  and  the  laboratory  temperature  in  degrees  Celsius.  The  correction 
equation  was  developed  from  information  provided  by  the  manufacturer.^^^) 

Corrected  Hook  Gage  Pressure  = 

(Observed  Hook  Gage  Pressure)(Correction  Factor)  . (7) 


Table  2.  Calibration  Summary 


Equipment 

Number  of  Calibration  Points 

0-4  inch  Ranee  4-6  inch  Range 

Transducer  No.8 

290 

0 

Transducer  No.5 

305 

16 

Transducer  No.4 

295 

0 

Transducer  No.  3 

252 

0 

Inclined  Manometer  No.  1 

305 

N/A 

Inclined  Manometer  No.2 

305 

N/A 

Meri-cal® 

1 38  points  in 

the  0-12  inch  range 
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Temperature  (®C) 


Figure  19:  Hook  Gage  Pressure  Correction  Factor  For  Indicating  Ruid  Density  Changes. 


ASCn  text  files  produced  by  the  data  acquisition  software  were  opened  and  read  by  the 
“Calibration  Digestion”  computer  program  written  by  Dr.  Guffey.  This  program 
determined  which  transducer  was  on  and  calculated  the  estimated  pressure  based  on 
previously  determined  calibration  coefficients.  The  resulting  values  were  recorded  to  a 
“Digested  Calibration  Data”  ASCII  file  containing  the  date,  time,  identification  code, 
transducer  number,  measured  voltage,  and  estimated  pressure  for  each  calibration  point. 

The  contents  of  the  calibration  data  file  were  imported  into  a  computer  spreadsheet 
program,  and  the  observed  pressure  values  from  the  hook  gage  and  other  instruments  were 
entered  adjacent  to  the  proper  identification  code.  Hook  gage  values  were  corrected  for 
temperature  using  the  relationship  developed  for  indicating  fluid  density  and  temperature. 
Next,  a  text  file  of  the  resulting  data  was  imported  into  a  commercial  statistics  program 
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(Data  Desk®  Professional,  Odesta  Corporation,  Nonhbrook  IL).  Scatter  plots  of  corrected 
hook  gage  pressures  versus  transducer  voltages  and  of  the  hook  gage  pressures  versus  the 
pressures  observed  with  other  instruments,  identified  errors  resulting  from  miskeyed  data, 
transposed  numbers  or  spurious  transducer  voltages.  The  sources  of  errors  were 
investigated  by  reviewing  the  original  data  collection  sheets.  When  the  source  of  error  was 
found  it  was  corrected;  however,  data  points  with  unexplained  error  were  removed  from 
the  data  set  for  analysis  since  they  would  exert  high  leverage  during  regression  analysis. 

Transducer  voltages  were  related  to  the  pressures  measured  with  the  hook  gage  using 
the  following  regression  equation: 

Hook  Gage  Pressure,  -  Cq  +  Ci (Voltage)  +  C2(Voltage2)  . (8) 

v.*iere;  Cq,  Ci,  and  C2  =  regression  coefficients  providing  the  best  fit  to 

the  pressure  range  of  interest 

Linear  regression  was  used  to  provide  calibration  curves  for  the  inclined  manometers 
and  the  Meri-cal®.  Following  is  the  general  form  of  the  model  used: 

Hook  Gage  Pressure  =  Cq  +  Ci (Pressure  Reading)  . (9) 

where:  Cq,  and  Ci  =  regression  coefficients  providing  the  best  fit  to 

the  pressure  range  of  interest 


The  pressure  transducers  apparently  experienced  smaU  daily  shifts  in  zero  and  gain. 
Figure  20  shows  an  apparent  zero  shift  for  one  set  of  data,  and  Figure  21  shows  an  apparent 
gain  shift  coupled  with  a  zero  shift.  As  a  result,  transducers  were  calibrated  for  each  set  of 
measurements  using  experimental  pressures  measured  with  an  inclined  manometer. 

Appendix  F  contains  the  calibration  curves  for  the  pressure  measuring  devices  used  in 
this  study. 
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Observed  F*ressure  (inches  w.g.) 


Figure  2 1 :  Illustration  of  Transducer  Zero  and  Gain  Shift 


Other  Measuring  Devices 

Barometric  pressure  was  measured  with  a  "National  Weather  Service  type"  barometer 
(Princo,  Model  453,  Southhampton  PA).  All  readings  were  corrected  for  temperature  and 
local  gravity  effects  using  equations  provided  by  the  manufacturer.(^5)  Corrected  observed 
values  from  the  laboratory  were  periodically  checked  against  corrected  values  provided  by 
the  local  U.S.  Weather  Bureau.  Readings  generally  differed  by  less  than  1  millimeter  of 
mercury. 
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Duct  temperatures  were  measured  with  a  digital  thermometer  (Fluke,  Model  52  K/J, 
Everett  WA).  After  a  15  minute  warm-up  period,  the  digital  readout  on  this  thermocouple 
was  compared  to  a  mercury  thermometer  (Ever  Ready  Thermometer  Co.,  NY)  that  was 
accurate  to  0.5°C  and  traceable  to  the  National  Bureau  of  Standards.  The  offset  for  the 
digital  thermometer  was  adjusted  if  necessary  so  that  the  digital  readout  was  within  0.2°C  of 
the  thermometer. 

Room  dry  bulb  and  wet  bulb  temperatures  were  measured  with  a  battery  powered 
psychrometer  (Cole-Parmer  Instrument  Company,  "Psychro-Dyne",  Model  3312-40, 
Chicago  IL).  The  dry  bulb  reading  for  both  thermometers  was  compared  t'^  the  mercury 
thermometer  mentioned  above.  A  discrepancy  of  less  than  2°C  was  considered  insignificant 
since  the  measurements  were  used  only  to  determine  the  contribution  of  water  vapor  to  the 
airstream  density  factor.  In  general,  measurements  were  within  0.5°C. 

Three  different  types  of  steel  measuring  tapes  were  used  to  measure  lengths  along  the 
ducts.  A  caliper  (manufactured  in  China,  manufacture  unknown)  with  a  range  up  to  9  inches 
(readable  to  0.(X)1  inch)  was  used  to  measure  the  diameters  of  the  ducts,  connectors,  and 
DOM  measuring  stations.  The  caliper  was  also  used  to  measure  the  gage  of  the  ducts  and 
connectors. 

An  electronic  tachometer  (Red  Lion  Controls,  "DITAK  5",  Model  DT-5,  York  PA)  was 
used  to  measure  the  fan  rotation  rate.  The  tachometer  determined  the  fan  rotation  rate  by 
counting  the  teeth  on  a  gear  mounted  directly  on  the  fan  drive  shaft.  It  was  readable  to  one 
rotation  per  mir  ute. 


Friction  Loss  Determination 

Friction  loss  equations  were  determined  for  4,  5,  6,  7,  and  10  inch  diameter  ducts  for 
velocities  ranging  from  1000-6000  feet  per  minute.  The  friction  loss  equations  for  5,  7  and 
10  inch  diameter  ducts  were  used  in  this  study. 


To  determine  friction  losses,  10  foot  lengths  (±3  inches)  of  spiral  wound  duct  were 
placed  on  the  ceiling  mounted  frame  described  earlier.  Appendix  H  contains  a  detailed 
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outline  of  the  procedure  used  to  collect  friction  data.  Figure  22  shows  a  generalized  setup 
for  all  friction  tests. 


Elevation  VP  Measurement  Station 

(DOM  tubing) 

0  30-40  D  1  2  3  4  ^56 


r 

*  i 

v  . ^  ' 

6'4" 

to  underside 
of  duct 

k 

j,.p“ . . .  E 

Figure  22:  Friction  Measurement  Setup 


After  assembly,  0.16  inch  (approximate)  diameter  holes  were  drilled  at  different 
locations  along  the  duct  (see  Figure  22  and  Table  3)  to  allow  centerline  static  pressure 
measurements.  The  fu’st  measuring  location  was  placed  approximately  30-40  duct 
diameters  downstream  from  the  duct  inlet  to  allow  the  friction  loss  per  foot  of  duct  to 
become  a  constant  value.  All  measurement  holes  were  placed  at  least  2  duct  diameters 
downstream  and  at  least  1  duct  diameter  upstream  from  duct  connectors. 

Short  sections  of  5  inch  diameter  duct  were  used  to  provide  the  desired  distance 
between  the  lateral  velocity  pressure  measurement  stations  and  the  junction  fitting.  As  a 
result,  two  connectors  were  required  between  the  DOM  tubing  and  the  junction.  Since  the 
friction  loss  per  foot  for  this  section  was  expected  to  differ  significantly  from  tests 
conducted  on  10  foot  sections,  this  short  section  of  duct  was  placed  after  the  DOM  tubing 
to  model  the  friction  losses  that  would  occur  during  junction  data  collection.  Table  4  shows 
the  number  of  connectors  per  test  section  for  each  setup. 
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Table  3:  Lengths  of  Test  Sections  for  Friction  Runs 


Length 

of  Test  Section  in  inches 

Diameter  (LL 

1:2 

hi 

1:1 

M 

M 

4  177 

321.5 

101.4 

422.9 

544.0 

121.1 

40.5 

5  208.5 

302.5 

99.8 

402.3 

513.3 

111.0 

48.3 

6  272.3 

98.5 

121.6 

220.1 

447.5 

227.5 

50.5 

7  258.5 

280.4 

101.6 

382.0 

497.3 

115.3 

46.5 

10  334.6 

X 

X 

X 

497.3 

257.5 

X 

Table  4: 

Number  of  Connectors  per  Test  Section  for  Friction  Runs 

Number  of  Connectors  per  Test  Section 

Diameter 

hi 

hi 

hi 

M 

li4 

M 

4 

3 

0 

3 

4 

1 

2 

5 

3 

0 

3 

4 

1 

2 

6 

0 

1 

1 

3 

2 

2 

7 

3 

0 

3 

4 

1 

1 

10 

X 

X 

2 

4? 

2? 

X 

A  duct  leak  test  was  accomplished  and  the  duct  alignment  was  verified  before  each  set 
of  data  was  collected.  The  fan  rotation  rate  was  adjusted  to  provide  velocities  in  the  range 
of  1000-6000  ft/min.  Velocity  pressures  from  the  three  10-point  traverses  were  measured 
with  a  pressure  transducer  and  the  resulting  voltages  were  recorded  to  a  computer  file. 
Centerline  velocity  pressures  were  also  measured  with  a  transducer  and  an  inclined 
manometer.  All  data  not  acquired  with  the  computer  were  recorded  to  a  friction  data 
collection  worksheet  (see  Appendix  I). 


The  static  pressure  at  duct  centerline,  the  airstream  temperature,  aiid  the  laboratory 
environmental  conditions  (wet  and  dry  bulb  temperatures  and  barometric  pressures)  were 
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recorded  and  used  to  calculate  the  airstream  density  factors.  The  centerline  static  pressure 
was  measured  with  the  transducer  and  a  redundant  measurement  was  taken  with  an  inclined 
manometer  (if  the  pressure  was  less  than  4  inches  w.g.)  or  with  the  Meri-cal®  (if  the 
pressure  was  greater  than  4  inches  w.g.).  Airstream  temperatures  were  measured  near  the 
DOM  tubing  with  the  digital  thermometer. 

All  static  pressure  differentials  were  measured  with  a  pressure  transducer  and  an  inclined 
manometer  using  Pitot  tubes  inserted  to  centerline  depths.  If  other  Pitot  tubes  were  between 
the  two  measuring  points,  they  were  raised  to  the  sidewall  of  the  duct  to  minimize  any  effect 
on  the  static  pressure  differential. 

The  "Digest  Workbench  Data"  computer  program  was  used  to  determine  the  average 
velocity  and  density  factor.  Occasionally,  static  pressures  exceeded  the  capacity  of  the 
transducer  or  were  in  error.  In  these  cases,  the  density  factor  and  predicted  average  velocity 
were  corrected  using  the  static  pressure  measured  with  the  inclined  manometer  or  the 
Meri-cal.® 

The  relationship  between  friction  loss  and  the  average  velocity  for  each  duct  was 
determined  by  conducting  linear  regression  analysis  on  the  log  of  the  differential  static 
pressures  measured  along  a  length  of  duct  and  the  log  of  the  corresponding  average 
airstream  velocity: 

log  (SPF/L)  =  Co  +  Cl  log  V  . (10) 

where:  SPF  =  differential  static  pressure  measured 

L  =  length  of  duct  between  measuring  location 
V  =  average  velocity  of  air  in  duct 
Cq  and  Cl  are  regression  coefficients 

The  differential  static  pressures  measured  with  the  inclined  manometer  were  used  in  these 
analyses  except  in  the  few  cases  where  the  manometer  value  was  clearly  inconectly 
recorded.  When  inclined  manometer  readings  were  deemed  incorrect,  the  differential  static 
pressures  predicted  by  the  transducers  v  ere  used  instead.  The  predicted  friction  per  length 
was  calculated  by  rearranging  Equation  10  as: 


SPF  =  l(  10^0  v^l) 
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where: 


SPF  =  differential  static  pressure  measured 
L  =  length  of  duct  between  measuring  location 
V  =  velocity  of  air  in  duct 
Q)  and  Cj  are  regression  coefficients 


Junction  Power  Loss  and  Modeling 

Power  Loss  Determination 

Many  of  the  methods  used  to  collect  junction  data  were  identical  to  methods  used  in 
collecting  ftiction  data  which  were  discussed  in  detail  earlier.  Only  the  methods  that  differ 
from  friction  are  discussed  in  detail  in  this  section.  Appendix  J  contains  a  complete  outline 
of  the  junction  data  collection  procedure. 

The  bilateral  junction  used  in  this  study  was  purchased  from  a  manufacturer  of 
ventilation  system  components  (United  McGill,  Stockton  CA).  The  junction  was 
constructed  from  18  gauge  galvanized  steel  that  was  rolled  into  separated  pieces  and  welded 
together  to  form  the  junction  (see  Figures  23  and  24).  The  laterals  were  welded 
approximately  midway  on  the  expansion  and  had  sharp  edges  at  the  point  of  connection. 
The  junction  was  designed  so  each  connection  was  slightly  smaller  than  the  nominal 
diameter  of  the  cormecting  duct;  this  allowed  the  ducts  to  slide  onto  the  junction  without 
requiring  a  separate  fitting.  The  inside  diameter  values  shown  in  Figure  24  are  the  average 
of  5  inside  diameter  readings  taken  with  the  caliper. 

Before  the  junction  was  used,  the  inside  surfaces  at  the  lateral  connections  were  ground 
with  a  grinding  wheel  attached  to  a  hand  drill  to  remove  all  metal  projections  in  the  path  of 
the  airstream.  Figure  25  is  a  photograph  of  the  inside  of  a  junction  before  grinding;  Figure 
26  shows  the  inside  of  a  junction  after  grinding. 


The  system  was  assembled  on  the  ceiling  mounted  frame  (see  Figures  27  and  28).  All 
ducts  and  measuring  stations  were  aligned  and  the  joints  were  sealed.  Figure  29  shows  the 
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system  components  and  the  general  layout  of  the  system  and  Figure  30  shows  all 
dimensions  for  the  system.  Figure  31  is  an  enlargement  of  the  junction  showing  the 
distance  from  the  measuring  station  to  the  junction  for  each  of  the  upstream  branches.  The 
dimensions  shown  in  Figure  31  and  the  distance  between  the  junction  and  the  measuring 
station  in  the  main  shown  in  Figure  30  were  used  to  calculate  the  lost  power  due  to  friction. 
Two  inches  were  added  to  the  distance  between  the  branch  measuring  stations  and  the 
junction  to  account  for  the  length  of  the  Pitot  probe  along  the  axjs  of  the  duct.  Similarly, 
two  inches  were  subtracted  from  the  distance  between  the  junction  and  tne  measunng 
station  in  the  main. 


Figure  23:  Photograph  of  Bilateral  Junction  Tested 
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Figure  25:  Inside  of  Junction  Before  Grinding 


Figure  26:  Inside  of  Junction  After  Grinding 


Branch "B” 


Figure  not  to  scale 


Figure  29:  System  Components  and  Layout 
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Figure  not  to  scale 
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Figure  30:  System  Dimensions 
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Data  were  collected  for  35  separate  runs  where  the  velocity  ratios  of  the  branches 
entering  the  junction  were  varied  from  0-3.7.  The  airflows  in  the  upstream  ducts  were  set 
using  the  impedance  method^'^^  to  achieve  the  desired  velocity  ratios.  Appendix  K  explains 
this  method  in  detail. 

Data  collection  and  reduction  procedures  for  each  branch  (and  main)  were  nearly 
identical  to  those  completed  for  friction  analyses.  The  only  difference  was  the  number  of 
differential  static  pressures  collected.  Detailed  methods  for  computer  data  acquisition  were 
discussed  earlier  in  the  "Computer  and  Data  Acquisition  and  Reduction"  section.  Appendix 
L  contains  a  copy  of  the  junction  data  collection  worksheet  used  to  record  data  not  collected 
with  the  data  acquisition  program. 

Final  data  reduction  was  performed  using  the  calibration  equations  for  the  transducers 
based  on  redundant  measurements  collected  with  the  inclined  manometer.  This  method  was 
used  because  of  transducer  zero  and  gain  shift  that  was  discussed  earlier.  The  priorities 
used  to  determine  which  static  pressure  measurements  should  be  used  in  data  analyses  were 
as  follows: 

1 .  First  priority  -  inclined  manometer  measurements,  unless  clearly  wrong 

2 .  Second  priority  -  transducer  values. 

3.  Meri-cal®  readings 

Once  the  data  was  reduced  to  its  final  form  all  relevant  data  was  merged  into  one 
computer  file.  Most  of  the  final  calculations  were  accomplished  on  a  computer  spreadsheet 
program.  Following  is  a  list  of  the  major  calculations  performed  on  the  spreadsheet  for 
each  junction  setup: 

1 .  Density  factor  correction 

2 .  Branch  velocity  ratios 

3 .  Average  velocity  for  each  branch 

4 .  Airflow  for  each  branch 

5 .  Static  pressure  at  each  measuring  station  (based  on  the  static  pressure  differential 
from  the  cross  section  in  duct  S) 

6.  "Mass"  airflow  in  the  main 

7 .  Velocity  in  the  main  based  on  the  "mass"  flow 
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8 .  Average  velocity  pressure  in  the  main  based  on  the"mass"  flow 

9.  Total  pressure  for  each  duct  based  on  the  "mass"  flow 

10.  Power  in  each  branch 

1 1 .  Total  lost  power 

1 2.  Lost  power  due  to  friction 

13.  Lost  power  of  the  junction 

The  static  pressure  measured  by  the  transducer  at  the  measuring  station  was  used  to 
calculate  the  airstream  density  factor.  The  density  factor  was  corrected  for  each  branch 
since  the  transducer  values  were  given  a  lower  priority  than  the  inclined  manometer  values. 
The  following  equation  was  used: 

cDFduci=DFduct  1.876^x71)  . ^  ^  ^  ^ 

where:  DF^uct  =  density  factor  computed  by  data  reduction  program 
cDFduct  =  corrected  density  factor 
Pbar  =  barometric  pressure 

SPused  =  static  pressure  at  the  measuring  station  used  (usually 
the  manometer  reading) 

SPTran  =  Static  pressure  at  the  measuring  station  sensed  by  the 
transducer 

The  factor  of  1.876  in  the  Equation  12  is  an  equivalence  factor  between  inches  of  water  and 
millimeters  of  mercury  (1  inch  water  =  1.876  mm  of  mercury).  This  conversion  factor  from 
the  ACGIH  Ventilation  ManualfO  was  later  found  to  be  incorrect  (1.865  is  the  correct 
value;(2)  fortunately  this  error  had  a  trivial  impact  on  the  data  analyses  since  the  difference 
produced  in  the  pressure  was  less  than  the  resolution  of  the  barometer  (0.1  mm  of 
mercury). 

The  corrected  average  velocity  for  each  duct  was  calculated  using  the  following 
formula: 


where:  cVavg  =  corrected  average  duct  velocity 

tVPavg  =  transducer  derived,  average  velocity  pressure 
calculated  by  the  data  reduction  program 


cDFduct  =  corrected  duct  density  factor 

The  static  pressure  differentials  measured  between  the  measuring  station  in  duct  S  and 
the  other  ducts  were  used  to  calculate  the  static  pressure  at  each  measuring  station.  This 
was  done  to  reduce  the  potential  errors  induced  by  slight  variations  in  the  fan  rotation  rate 
over  time.  The  differential  pressures  were  used  in  the  following  equations  to  give  the  static 
pressure  for  each  cross-section: 


calSPs  =SPs  (14) 

calSPn,  =  calSPj  +  ASPj.n,  (15) 

calSPa  =calSPs-ASPs.a  (16) 

calSPb  =  calSPs  -  ASPs.b  (17) 

where:  calSP  =  calculated  static  pressure  at  the  measuring  station 

ASPs.,  =  s*^6c  pressure  differential  measured  between  duct  S 
and  duct  i 


Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in 
ducts  A,  B,  S  and  M 

Similarly, the  sum  of  the  flows  in  the  branches  was  used  as  the  "mass"  flow  in  the 
main.  Power  loss  calculations  were  based  on  the  "mass"  flow  rates.  The  measured  values 
of  the  flow  in  the  main  were  compared  with  the  "mass"  flows  to  identify  any  errors. 

The  total  lost  power  was  calculated  from  the  sum  of  the  powers  of  the  upstream 
branches  and  the  power  in  the  main  using  Equation  2.  Using  the  nomenclature  in  Equations 
14-16  above.  Equation  2  becomes: 


43 


LP  =  Qa(caISPa  +  tVPavg  a)  +  ObCcalSP^  +  tVPavg  b)  + 

Q,(calSPs  +  VPavg  s)  -  Qm(calSPm  +  tVPav. . ^  . (18) 

where:  LP  =  total  lost  power 

Q  =  flow  in  the  duct 

calSP  =  calculated  static  pressure  in  the  duct 

tVPavg  =  average  velocity  pressure  calculated  by  the  data 
reduction  program 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in 
ducts  A,  B,  S  and  M 

As  discussed  in  the  Introduction,  the  total  lost  power  has  two  components,  the  power 
dissipated  due  to  the  airflows  mixing  and  the  power  dissipated  because  of  friction  between 
the  measuring  stations  and  the  junction.  The  estimated  static  pressure  loss  due  to  friction 
was  calculated  using  the  loss  equations  developed  in  the  fnction  study,  the  average  velocity 
in  the  duct,  and  the  distance  from  the  measuring  station  to  the  junction.  The  lost  power 
resulting  from  friction  in  all  of  the  ducts  was  estimated  using  Equation  3. 

Finally  the  lost  power  due  to  the  junction  was  found  by  subtracting  the  lost  power 
attributable  to  friction  from  the  total  lost  power  using  Equation  5.  Using  the  nomenclature 
from  Equations  14-18  above,  Equation  5  transforms  to  the  following: 

LPJ  =  [(3a(calSPa  +  tVPavg  a)  +  Qb(calSPb  +  tVPavg  b)  . (19) 

+  Qs(calSPs  +  VPavg  s)  -  Qm(calSPm  +  tVPavg  m)] 

-  [Qm(ASPFj.n,)  Qa(ASPFj.a)  +  Qb(ASPFj.b)  +  Qs(ASPFj.s)] 

where:  LPJ  =  {X)wer  lost  due  to  the  junction 

ASPFj.j  =  predicted  change  in  static  pressure  (due  to  friction) 
between  the  junction  and  the  duct  measuring  station 

Q  =  airflow  in  the  duct 

calSP  =  calculated  static  pressure  in  the  duct 

tVPavg  =  average  velocity  pressure  calculated  by  the  "Digest 
Workbench  Data  program 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in 
ducts  A,  B,  S  and  M 
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Junction  Modelling 

Multiple  regression  analysis  was  accomplished  using  the  kinetic  powers  in  the  upstream 
branches  as  the  independent  variables  and  the  lost  power  due  to  the  junction  (LPJ)  as  the 
dependent  variable  (see  Equation  b)/"*)  Unlike  the  approach  used  by  Guffey  and  Fraser,  the 
kinetic  power  of  the  main  was  not  used  as  an  independent  variable  since  it  is  not 
independent  of  the  kinetic  powers  in  the  upstream  ducts. 

The  wide  range  of  velocity  ratios  tested  was  important  in  evaluating  the  robustness  of 
the  kinetic  power  model.  Some  ratios  were  extreme  enough  that  a  poor  fit  to  the  kinetic 
power  model  was  expected.  Under  extreme  conditions  the  flow  decelerates  from  a  high 
velocity  to  a  sharply  lower  velocity  in  the  main.  In  these  cases,  the  power  losses  may  be 
more  a  response  to  a  sudden  expansion  than  a  function  of  mixing  phenomena.  To  evaluate 
model  performance  under  more  realistic  conditions,  subsets  of  data  were  analyzed  using 
narrower  velocity  ratio  ranges. 

Regression  diagnostics  were  executed  to  check  multiple  regression  assumptions  and  to 
identify  potential  outliers  in  the  data.  Table  5  provides  a  summary  of  the  regression 
diagnostics  used  and  their  purposes. 


Table  5:  Regression  Diagnostics  Summary 


Diagnostic  Tool 
Partial  Regression  Plot 

Leverage  Histogram 

Externally  Studentized  Residuals 

Cook's  Distance  and  DFFITS 


Tests  for 

Linearity  between  each  independent  variable  and 
the  dependent  variable. 

How  extreme  the  independent  variables  are  in  the 
individual  test  cases  . 

How  extreme  the  dependent  variable  is  in  the 
individual  test  cases. 

The  overall  influence  of  the  individual  test  cases  on 
the  regression. 


CHAPTER  III 

DISCUSSION  OF  RESULTS 


Friction  and  Flow  Estimation 

Pressure  differentials  measured  across  ventilation  system  components  generally  include 
substantial  contributions  due  to  friction  losses  since  experimental  measurements  must  be 
taken  far  downstream  of  the  component  to  avoid  disturbances  induced  by  the  component 
being  evaluated.  This  is  particularly  true  for  junctions.  To  determine  the  contribution  due  to 
a  junction  alone,  it  was  necessary  to  subtract  the  estimated  contribution  of  friction. 
Similarly,  the  methodologies  most  widely  used  to  design  new  ventilation  systems 
(ASHRAE  and  ACGIH)  compute  dynamic  and  friction  losses  separately. 

An  accurate  estimation  of  airflow  in  each  duct  is  critical  for  the  total  power  loss 
calculations  and  for  estimating  the  power  lost  due  to  friction.  In  addition,  deviations 
between  the  measured  downstream  flow  and  the  sum  of  the  measured  upstream  flows 
would  produce  a  spurious  change  in  power  losses  calculated  using  Equation  2.  For  that 
reason,  the  sum  of  flows  in  the  ducts  upstream  of  the  junction  (the  "mass”  flow)  was  used 
as  the  flow  in  duct  M  (the  downstream  duct)  for  all  analyses. 

Friction  Estimation 

Estimating  the  contribution  from  friction  would  seem  straight  forward:  experimentally 
determine  the  friction  loss  per  length,  then  use  that  value  to  predict  the  static  pressure  loss 
in  the  experimental  setup  resulting  from  friction.  Unfortunately,  the  friction  loss  per  length 
varies  with  the  distance  downstream  of  a  disturbance.  For  example,  in  a  duct  containing  a 
bellmouth  hood  op)ening,  the  loss  per  foot  of  duct  does  not  become  a  constant  value  until  a 
distance  approximately  15  duct  diameters  downstream  of  where  the  hood  connects  to  the 
duct.(*^)  Therefore,  the  dynamic  loss  of  a  system  component  is  confounded  with  the  effects 
of  disturbed  flow  on  friction  losses. 

Error  is  introduced  by  applying  empirical  friction  loss  equations  based  on 
measurements  far  downstream  of  a  disturbance  to  a  section  of  duct  immediately 
downstream  of  a  disturbance.  However,  there  is  no  model  currently  available  to  estimate 
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friction  losses  in  these  circumstances.  As  a  result,  the  fiiction  losses  in  the  ducts  upstream 
and  downstream  of  the  junction  were  estimated  from  empirical  equations  that  were  based 
on  data  collected  far  downstream  from  the  inlet  of  a  long,  straight  duct  (see  Table  6  and 
Figure  32).  Appendix  M  contains  the  results  from  the  regression  analyses  conducted  to 
determine  these  relationships. 

Since  differential  pressures  were  measured  at  several  cross-sections  along  the  length  of 
a  single  duct,  it  was  possible  to  compare  the  friction  loss  per  foot  in  the  most  upstream  test 
section  (section  l<->2  in  Figure  22)  with  the  most  downstream  test  section  (section  3<-^4) . 
This  helped  establish  that  the  friction  per  unit  length  was  the  same  in  both  sections.  As  a 
result,  the  data  from  the  entire  length  (section  l<-^4)  was  used  to  calculate  the  friction  loss 
equations,  thus  minimizing  the  effects  of  measurement  error. 

Duct  connectors  are  responsible  for  a  p)ortion  of  the  overall  friction  loss  along  any 
length  of  duct  and  their  loss  per  unit  length  is  expected  to  be  higher  than  that  of  the  duct.  In 
this  study,  the  portion  of  total  friction  losses  attributable  to  the  connectors  was  estimated  by 
comparing  friction  losses  (on  a  per  foot  basis)  between  a  section  with  no  connectors  and  a 
section  containing  one  or  more  connectors.  Figure  33  shows  the  effect  of  varying  the 
number  of  connectors  in  a  100  foot  run  of  7  inch  diameter  duct.  The  effect  of  connectors 
was  similar  for  other  duct  diameters. 

The  short  section  of  duct  that  was  required  between  the  measuring  stations  and  the 
junction  for  the  5  inch  diameter  laterals  was  placed  after  the  DOM  measuring  station  during 
friction  data  collection  (test  section  5<->6  in  Figure  22)  to  model  the  increase  in  friction  loss 
per  unit  length  due  to  the  connector  effect  discussed  above.  Likewise,  the  DOM  connection 
to  the  junction  for  the  7  inch  diameter  duct  was  also  modeled  to  determine  any  increase  in 
friction  loss  per  unit  length  for  this  arrangement.  For  the  5  inch  diameter  duct  the  friction 
losses,  when  normalized  for  length,  were  approximately  1 .7  times  higher  in  test  section 
5<-^6  than  the  longest  test  section  ( 1  <-^4).  For  the  7  inch  diameter,  the  friction  losses  were 
only  1.3  times  higher  for  in  test  section  5<->6  than  section  l<->4.. 


Table  6:  Empirical  Friction  Loss  Equations 


Figure  32:.  Power  Law  Regression  F  o  Observed  Friction  Loss/100  feet  (Section  l<->4  ) 
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Figure  33:  Effect  of  the  Number  of  Connectors  on  Friction  Loss 


Accurate  empirical  models  for  friction  losses  were  very  important  since  the  lost  power 
due  to  the  junction  was  estimated  by  subtracting  the  lost  power  due  to  friction  from  the  total 
lost  power.  Analyses  revealed  the  power  lost  due  to  friction  accounted  for  a  significant 
portion  (an  average  of  58  percent)  of  the  total  power  lost  (see  Figure  34).  The  proponion  of 
power  losses  due  to  friction  were  highest  when  the  duct  airstream  velocities  were  similar. 
This  was  true  because  the  losses  due  to  the  junction  were  lowest  for  these  conditions.  Most 
of  the  power  losses  due  to  friction  occurred  in  the  relatively  long  section  of  duct 
downstream  of  the  junction  fitting.  Duct  M  (see  Figures  29  and  35).  The  contributions  due 
to  the  upstream  ducts  varied  as  a  function  of  the  velocities  through  them,  and  were 
relatively  small  in  all  cases. 
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Figure  34;  Percent  of  Total  Power  Loss  Resulting  from  Friction 
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Figure  35:  Contribution  of  Different  Ducts  to  the  Friction  Power  Losses 


Row  Estimation 

The  flows  in  the  downstream  main  (QnO  were  estimated  from  measured  velocity 
pressures  in  the  main  (Qmeas)  and  by  summing  the  upstream  airflows  (Qmass)-  In  twenty-six 
of  32  cases,  the  two  estimates  varied  by  less  than  1%,  but  8  varied  by  1-2%  and  one  varied 
by  3.1%  (see  Figure  36).  Figure  37  shows  that  the  percent  error  was  normally  distributed 
except  for  case  number  2  which  may  represent  a  measurement  error. 

Since  the  dependent  variable  (lost  power  due  to  the  junction  )  is  a  function  of  the 
difference  in  energy  rates  upstream  and  downstream  of  the  junction,  Qmass  was  used  for  (^m 
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for  power  computations  (see  below).  Otherwise  a  spurious  power  difference  would  exist 
attributable  solely  to  differences  between  Qmass  and  Qmeas- 


Figure  36:  Percentage  of  Error  Between  Measured  and  "Mass"  Flows  In 

Duct  M 
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Junction  Modelling 

“Observed”  values  of  lost  power  in  the  junction  (LPJ)  were  fit  to  the  linear  kinetic 
power  model,  which  was  modified  by  the  addition  of  a  regression  constant  (Co)  for 
purposes  of  statistical  regression: 

LPJ  =  QaVPaKa  +  QbVPfaKb  +  QsVPgKs  +  Co  . (20) 

Table  7  presents  the  results  of  several  multiple  regression  runs.  First  the  entire  range  of 
cases  were  analyzed,  then  various  subgroups  of  the  entire  data  set  were  selected  based  on 
velocity  ratio  ranges  and  were  analyzed  separately.  To  represent  physical  reality,  the  model 
should  satisfy  four  criteria  for  each  subset  of  data:  1)  it  should  fit  observed  data  well,  2)  at 
zero  flow,  it  should  predict  zero  losses  (Co=0),  3)  the  coefficient  values  Kg  and  Kb  should 
be  approximately  equal,  and,  4)  it  should  satisfy  diagnostic  criteria  for  applicability  of 
multiple  regression  analysis.  To  be  useful,  the  model  need  not  fit  the  entire  range  well  as 
long  as  its  fits  the  range  of  practical  application  well. 

Goodness  of  Fit 

Although  the  model  fit  the  full  range  of  data  moderately  well  (Adj  =  0.896),  its  fit 
improved  dramatically  (R2=0.977)  when  the  data  were  narrowed  to  the  velocity  ratios  that 
might  be  expected  of  high  velocity  ventilation  systems  (0.75<Vratio<l-25).  The  model  fit 
all  subsets  of  the  data  at  least  moderately  well  (R2>0.86).  The  values  of  R^  were 
determined  with  the  regression  constant  included,  but  the  values  of  Kg,  Kb,  and  Ks  in  Table 
7  were  determined  with  the  regression  constant  excluded  since  the  value  of  LPJ  must  be 
zero  at  zero  flow. 

Predicted  Value  at  Zero  Flow 

Since  no  losses  are  possible  at  zero  flow,  the  model  should  predict  a  zero  value  of  LPJ 
when  flows  are  zero.  In  fact,  the  regression  constant  for  the  "all  data"  category  was 
approximately  5%  of  the  mean  value  of  LPJ  while  that  for  the  "0.5<Vratio  <2.0  was  less 
than  1%  of  the  mean  value  of  LPJ  ;  these  v.'ere  excellent  findings.  On  the  other  hand,  the 
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regression  constant  for  the  ”0.75<Vratio<l-25"  was  approximately  24%  of  the  mean  value 
of  LPJ;  however,  this  finding  is  not  surprising  because  of  the  reduced  number  of 
observations. 


Table  7:  Regression  Results  for  Junction  Modelling 


Selector 

K 

Kai 

m 

Ksi 

Adi  R2** 

All  Data 

35 

0.613 

0.489 

0.188 

89.6% 

63.7 

0.5<Vratio‘<2.0 

19 

0.604 

0.374 

0.216 

90.9% 

36.0 

0.75<Vratio<1.25 

9 

0.494 

0.361 

0.254 

97.7% 

23.0 

Va=Vb 

16 

X 

1.10 

0.215 

95.7% 

43.9 

Va=Vb 

16 

1.06 

X 

0.215 

92.7% 

57.5 

Va^tVb 

19 

0.642 

0.579 

0.119 

94.1% 

102.6 

Va?tVs 

28 

0.635 

0.452 

0.199 

93.2% 

56.8 

Vb^Vs 

26 

0.597 

0.512 

0.197 

89.9% 

65.2 

Va^Vs5tVb 

9 

0.646 

0.569 

0.133 

86.4% 

61.4 

Vratio  =  Velocity  ratios 

iwilh 

no  regression  constant 

"with  regression  constant 

included 

Ideally,  the  value  of  Ka  would  be  exactly  equal  to  Kb  since  ducts  A  and  B  are  physically 
identical.  The  difference  between  Ka  and  Kb  may  occur  because  the  flows  in  the  upstream 
branches  were  not  set  randomly.  For  much  of  the  data,  the  velocity  in  two  of  the  three 
upstream  branches  were  set  equal  to  each  other  while  the  flow  in  the  third  was  varied. 

Thus,  Qa  and  VPa  were  completely  dependent  of  Qb  and  VPb;  the  intercorrelation  of 
independent  variables  is  shown  in  Table  8.  Intercorrelation  of  independent  variables 
produces  instability  and  distortion  of  the  regression  coefficients  even  when  the  coefficients 
accurately  predict  the  losses^'®-^^).  To  remove  the  intercorrelation,  a  subset  of  data 
(Va^tVb)  was  analyzed  to  produce  values  of  coefficients  for  the  kinetic  power  in  duct  A  and 
duct  B  that  were  much  closer  (see  Table  7).  Another  subset  (Va=Vb)  was  analyzed  by 
removing  the  kinetic  power  in  duct  A  as  a  variable  and  then  the  kinetic  power  in  duct  B  as 
a  variable.  When  this  was  done  the  values  of  the  coefficients  were  also  much  closer  (see 


Table  7).  This  further  supports  the  likelihood  that  the  differences  in  the  coefficients  for 
ducts  A  and  B  are  an  artifact  of  the  multiple  regression  fitting. 


Table  8:  Summary  Analyses  for  Collinearity  Between  Independent  Variables 


First  Variable 

Second  Variable 

R2 

Kinetic  Power  in  duct  S 

Kinetic  Power  in  duct  A 

14.6% 

Kinetic  Power  in  duct  B 

Kinetic  Power  in  duct  A 

27.1% 

Kinetic  Power  in  duct  B 

Kinetic  Power  in  duct  S 

36.8% 

Regression  Diagnostics 

As  will  be  shown  next,  except  for  the  for  the  slight  collinearity  of  the  independent 
variables  discussed  above  the  basic  multiple  regression  assumptions  were  met: 

1 .  Independent  variables  are  truly  independent  of  each  other.  As  discussed  earlier,  a 
check  for  collinearity  indicated  there  was  a  slight  relationship  between  the 
independent  variables.  See  Table  8. 

2 .  There  is  a  linear  relationship  between  each  independent  variable  and  the  dependent 
variable. 

3 .  No  value  of  a  independent  or  dependent  variable  exerts  undue  influence  on  the 
regression  analysis. 

The  linear  relationship  between  the  lost  power  in  the  junction  and  the  kinetic  power  in 
each  branch  was  verified  by  creating  partial  regression  plots  which  show  the  relationship 
between  the  dejjendent  variable  and  an  independent  variable  with  the  effects  from  the  other 
independent  variables  removed.  The  partial  regression  plots  showed  that  a  linear 
relationship  existed  between  the  lost  power  in  the  Junction  and  each  of  the  kinetic  powers  in 
the  branches 
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The  influence  of  the  particular  values  of  independent  and  dependent  variables  was 
checked  using  leverage  values  and  externally  studentized  residuals,  respectively.  A 
histogram  of  leverages  indicated  that  case  numbers  2  ,4,  and  6  contained  extreme  kinetic 
power  values.  A  histogram  and  normal  probability  plot  of  externally  studentized  residuals 
indicated  that  case  number  5  had  an  extreme  junction  power  loss  value  as  compared  with 
the  other  case. 

Histograms  of  Cook's  distances  and  DFFITS  were  used  as  diagnostic  tools  to  assess 
the  overall  influence  of  any  one  case  on  the  regression  results.  These  tests  indicated  that 
case  numbers  5  and  6  may  exert  undue  influence  on  the  regression  results. 

The  influence  of  these  possible  outliers  was  investigated  by  removing  each  case  number 
identified  from  the  regression  analysis.  Removing  case  numbers  5  and  6  from  the 
regression  improved  the  regression  without  having  a  substantial  impact  on  the  coefficients. 
Table  9  provides  a  summary  of  the  results  of  all  analyses. 

The  data  for  the  cases  identified  as  possible  outliers  were  reviewed  to  see  if  removing 
these  cases  from  the  regression  was  justified.  Histograms  of  the  kinetic  powers  identified 
that  case  numbers  5  and  6  were  the  only  cases  that  had  a  combination  of  both  the  highest 
and  lowest  kinetic  powers  of  all  the  case  numbers.  Case  number  5  had  the  highest  kinetic 
power  in  duct  B  of  all  cases  while  also  having  extremely  low  kinetic  powers  in  ducts  A  and 
S.  Similarly,  case  number  6  had  the  highest  kinetic  power  in  duct  S,  while  the  kinetic 
powers  in  ducts  A  and  B  were  among  the  lowest.  Since  case  numbers  5  and  6  were 
extreme  because  of  experimental  design  and  not  because  of  error,  there  was  no  justification 
for  removing  them  when  reporting  the  final  results,  especially  since  removing  them  had 
negligible  effect  on  the  values  predicted  by  the  model. 
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Table  9;  Summary  of  Regression  Results  from  After  Removing  Possible  Outlying  Data 


Case  .\umber(s) 

Removed 

Regression  Coefficients 

Ka  Kb  Ks 

Adj  R2 

se 

None 

0.613 

0.489 

0.188 

89.6 

63.7 

2 

0.583 

0.503 

0.188 

89.8 

63.3 

4 

0.602 

0.475 

0.193 

89.9 

62.8 

5 

0.628 

0.446 

0.197 

92.0 

55.8 

6 

0.617 

0.511 

0.169 

90.6 

60.6 

2, 4, 5,6 

0.584 

0.455 

0.192 

94.0 

48.3 

5,6 

0.692 

0.467 

0.182 

93.7 

53.89 

Influence  of  Branch  Velocity  Ratios  on  Junction  Lost  Powers 

Figures  38—43  shoves  the  relationship  between  the  normalized  lost  power  in  the 
junction  and  the  vai^cus  velocity  ratios.  For  graphical  presentation,  the  lost  power  in  the 
Junction  for  each  case  number  was  normalized  by  the  kinetic  power  in  duct  M,  thus 
reducing  each  junction  loss  to  an  equivalent  velocity  pressure  coefficient  in  the  downstream 
duct.  As  expected,  the  normalized  lost  power  due  to  the  junction  (LPJ/QmVPm)  was  lowest 
when  the  velocities  of  all  three  incoming  branches  were  approximately  equal  (see  Figure 
38).  Large  velocity  differences  at  the  Junction  were  expected  to  increase  separation  from  the 
walls  due  to  momentum  or  sudden  expansions,  thus  increasing  power  losses. 

Velocity  in  the  Side  Branches  With  Respect  to  the  Center  Branch 

Figures  38  shows  the  relationship  between  the  normalized  lost  power  in  the  junction 
and  velocity  ratio,  VbA's.  Ttie  apparent  outliers  in  Figure  38  result  from  extreme  values  of 
the  Va  compared  to  V^,  and  Vs.  Figure  39  shows  the  normalized  LPJ  with  the  effects  of  the 
VaA'b  and  Va/V s  ratios  removed  by  plotting  only  the  ca.ses  where  the  velocities  in  ducts  A 
and  B  were  equal. 
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LPJ/Qm  VPm 

0. 

0. 
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Figure  38:  Influence  of  Velocity  Ratio  (Vb/Vj)  on  Junction  Lost  Power 

Figure  39  shows  that  power  losses  were  smallest  when  the  velocities  in  all  three  ducts 
entering  the  junction  were  approximately  equal.  Interestingly,  this  figure  also  shows  that 
power  losses  increased  slightly  when  the  decreased  with  respect  to  Vs,  but  increased 
rapidly  as  Vb  exceeded  Vj.  This  differential  effect  was  expected  since  the  momentum  of 
flows  in  the  lateral  ducts  would  increase  separation,  while  the  momentum  of  flows  in  duct 
S  would  tend  to  reduce  separation  of  flows  A  and  B.  The  relationship  between  the  power 
lost  in  the  junction  and  the  ratio  Vg/Vj  was  similar  (see  Figures  40  and  41). 


0  0.5  1  1.5  2  2.5 

Vb/Vs 


Figure  40:  Influence  of  Velocity  Ratio  (Va/Vj)  on  Junction  Lost  Power 


Figure  42  shows  that  the  lost  power  in  the  junction  varies  with  respect  to  the  velocity 
ratio  of  the  laterals  (VaA^b).  The  independent  effects  of  this  ratio  are  not  clear  due  to  the 
confounding  effects  of  the  velocity  ratios  and  Vi/Vs-  Figure  43  shows  that  the  power 
losses  rose  rapidly  when  Vb  increases  with  respect  to  Vg  and  Vs.  The  two  points  directly 


above  the  Vb/Vg  =  1.0  in  Figure  43  had  nearly  the  same  velocity  ratios  but  different  levels 
of  total  air  flow.  Ideally,  these  two  points  would  have  approximately  the  same  normalized 
lost  powers.  The  difference  is  likely  due,  in  part,  to  a  slight  difference  in  velocity  ratios.  All 
velocity  ratios  were  rounded  to  nearest  0.1;  therefore,  the  true  velocity  ratios  could  differ 
by  nearly  10%  and  still  appear  equal  on  the  graphs.  Some  of  the  deviation  also  likely 
results  from  differing  errors  in  estimating  friction  and  other  measurement  errors. 


LPJ/Qm  VPm 


I 


Figure  42:  Influence  of  Velocity  Ratio  (Vt/Va)  on  Junction  Lost  Power 


a 


LPJ/Qm  VPm 


WVa 


Figure  43:  Influence  of  Velocity  Ratio  (Vb/Va)  on  Junction  Lost  Power  Vg  =  Vb 


Estimation  of  "True"  Coefficients  for  the  Junction  Tested 

As  discussed  in  the  junction  modelling  and  regression  diagnostics  sections,  the 
experimental  design  produced  a  degree  of  collinearity  of  independent  variables,  which  lead 
to  non-equal  values  for  the  regression  coefficients  Kg  and  Kb.  Since  it  is  misleading  to 
assign  different  values  to  two  identical  branches,  the  "true"  (equal)  values  of  Kg  and  Kb 
w’ere  each  estimated  by  averaging  Kg  and  Kb- 
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WVs 

Figure  44:  Comparison  of  Kinetic  Power  Coefficients  Provided  by  Regression  and  an 

Average  of  Ka  and  Kb 

Figure  44  shows  a  close  correspondence  between  the  lost  powers  predicted  by  the 
regression  coefficients  and  those  predicted  using  the  averaged  values  of  Ka  and  Kb. 
However,  Figure  45  shows  there  is  some  error  when  comparing  either  set  of  predicted 
values  to  the  measured  values.  It  should  be  reemphasized  in  Figure  45  that  Va  was  not 
exactly  equal  to  Vb  in  any  case,  and  that  the  ratio  VbA'a  varied  somewhat  from  case  to  case 
due  to  the  rounding  of  velocity  ratios  as  discussed  earlier.  Since  the  predicted  and  observed 
values  of  LPJ  are  functions  of  VbA^a.  even  these  slight  variations  prevent  the  values  from 
following  the  smooth  curve  that  might  normally  be  expected. 
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Comparison  of  Power  Losses  From  the  Test  Bilateral  Junction  to  Losses 
from  an  Equivalent  Double  Single-Lateral  Junction  Setup 

The  second  goal  of  this  research  was  to  determine  if  the  power  losses  in  the  bilateral 
junction  were  significantly  higher  than  the  power  losses  predicted  for  an  equivalent  set  of 
two  single-lateral  junctions.  Figure  46  illustrates  a  hypothetical  layout  that  is  equivalent  to 
the  bilateral  junction  tested. 


Figure  46:  Hypothetical  Setup  with  Two  Single-Lateral  Junctions  Equivalent  to  the 

Bilateral  Junction  Tested 


The  power  losses  for  two  single-lateral  junctions  were  predicted  applying  kinetic 
powers  measured  in  the  branches  of  the  bilateral  junction  to  equations  develop)ed  by  Guffey 
and  Fraser^'*)  for  single-lateral  junctions.  The  following  equation  for  the  double  single¬ 
single-lateral  junction  setup  was  derived  by  combining  the  predicted  powers  for  two 
separate  junctions  into  one  expression: 

LPJpred  =  0.045  KPs  +  0.351  KPg  +  0.062  KPx  -i-  0.321  KPb  -t-0.017  KP^  . (21 ) 
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Use  of  Equation  20  assumes  "ideal"  conditions  of  no  friction  loss  between  the  junctions, 
no  interaction  effects  between  the  junctions,  and  no  additional  pxjwer  losses  due  to  longer 
laterals. 

The  measured  power  losses  for  bilateral  junctions  were  an  average  45%  higher  than  the 
power  losses  predicted  for  the  two  single-lateral  junctions.  The  range  of  the  excess  power 
losses  for  the  bilateral  junctions  across  all  conditions  was  14-92%.  The  excess  power 
losses  were  at  a  minimum  when  the  velocity  ratios  were  nearly  1 .0.  The  largest  difference 
(92%)  occurred  when  one  branch  of  the  bilateral  junction  was  completely  blocked  so  the 
junction  was  forced  to  behave  as  a  single-lateral  junction.  Table  10  is  a  summary  of 
analyses  demonstrating  that  the  excess  power  losses  were  highest  for  the  extreme  velocity 
ratios  tested  and  that  the  excess  power  losses  decreased  as  the  velocity  ratios  approached 
1 .0.  Figures  47  and  48  compare  the  normalized  measured  power  losses  in  the  bilateral 
junction  with  those  predicted  for  two  single-lateral  junctions. 


Table  10:  Summary  Statistics  for  the  Percentage  of  Excess  Power  Losses  in  the  Bilateral 
Junction  as  Compared  with  Predicted  Power  Losses  for  Two  Single  Lateral  Junctions 


Selector 

Avgragg  (%) 

Range  (%) 

All  data 

45.0 

14-92 

0.5<  V  ratio '^2.0 

32.3 

14-59 

0.75<Vratio<l-25 

27.1 

14-38 

Figure  47:  Comparison  Measured  Bilateral  and  Two  Single-lateral  Junction  Power  Losses 

for  Velocity  Ratio  (Vt/Vs)  Where  Vg  =  Vh 


LPJ/Qm  VPm 


Figure  48:  Comparison  of  Measured  Bilateral  and  Predicted  Two  Single-lateral  Junction 
Power  Losses  for  Velocity  Ratio  (Vb/Vg)  Where  Va  = 


$ii^..if!C^.nce  of  the  Higher  Power  Losses  For  the  Bilateral  Junction  Compared  to  the 
Single-lateral  Junction  Setup 
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Although  the  power  losses  appear  substantially  higher  than  those  predicted  for  two 
single-lateral  junctions,  it  is  important  to  place  these  increases  in  the  perspective  of  overall 
system  pressure  and  power  requirements.  To  provide  a  basis  of  comparison,  an  additional 
parameter  will  be  developed.  From  Equation  5  (assuming  no  friction  losses): 

QmTPm  =  QaTPs  +  QbTPb  +  QsTPs-LPJ  . (22) 

Where:  LPJ  =  lost  power  due  to  the  junction 

TP  =  total  pressure  in  a  duct 
Q  =  flow  in  a  duct 

Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in  ducts 
A,  B,  S  and  M 


Rearranging  equation  21  becomes: 


.(23) 


The  units  of  LPJ/Qm  are  pressure,  so  it  is  possible  to  define  TPJioss  as: 
LPJ 


Qm 


.(24) 


Where:  TPJiqsj  =  total  pressure  lost  due  to  the  junction 


If  both  sides  of  Equation  23  are  divided  by  the  downstream  velocity  pressure  (VPm),  then 
there  exists  an  equivalent  downstream  velocity  pressure  coefficient  such  that: 


py  _  2^PJ|qss 
^ "  VPm 


(25) 


Where: 


VP  =  velocity  pressure  in  a  duct 

Fj  =  equivalent  downstream  pressure  coefficient  for  a 
junction 
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Subscripts  a,  b,  s,  and  m  represent  reference  cross-sections  in  ducts 
A,  B,  S  and  M 


By  normalizing  lost  power  with  downstream  kinetic  power  to  produce  Fj,  the  power  loss 
due  to  the  mixing  in  the  junction  can  be  compared  to  loss  coefficients  for  elbows  and  other 
components  in  the  downstream  duct. 

As  shown  in  Figu.e  47,  the  value  of  Fj  for  a  bilateral  junction  ranges  from 
approximately  0.30  to  0.85  as  compared  to  a  range  of  about  0.2  to  0.5  for  two  single- 
lateral  junctions.  The  maximum  difference  of  approximately  0.35  is  less  than  that  of  two 
90°  elbows  (see  Figure  49).  Also  shown  are  equivalent  lengths  of  10  inch  diameter  ducts 
(size  of  main  leaving  the  bilateral  and  second  of  the  two  single  lateral  junctions).  Table  1 1 
presents  a  comparison  of  Fj  to  other  typical  loss  coefficients. 


Table  11:  Excess  Bilateral  Junction  Losses  as  Equivalent  Lengths  of  10  inch  Diameter  Duct 
and  90°  elbows  for  Velocities  in  the  Range  of  1000-4000  feet  per  minute 


4FJ* 

Number  of  Diameters 

Number  of  Elbows  ** 

0.30 

16-18 

1.6 

0.20 

11-12 

1.1 

0.15 

7-8 

0.79 

0.10 

4-5 

0.52 

0.05 

2-4 

0.26 

»p  _  LPJobserved  bilateral  ~  LPJpredicted  2-l3lerals 

“  QmVPm 


AFj  is  the  loss  coefficient  for  the  excess  power  losses  in  the  bilateral  junction  as 
compared  to  two  single-lateral  junctions. 


**  Elbows  are  5-section  with  a  radius  of  curvature  equal  to  twice  the  diameter. 


Finally,  it  should  be  noted  that  the  predicted  excess  pressure  losses  due  to  using  the 
bilateral  Junction  instead  of  two  single-lateral  junctions  would  be  very  small  in  comparison 
to  overali  losses  in  typical  ventilation  systems.  Even  for  the  simple  system  constructed  to 
test  the  bilateral  junction,  the  greatest  increase  in  total  pressure  was  only  3.5%  (Table  12), 
and  that  was  for  impractically  extreme  velocity  ratios. 

Table  12:  Percentage  of  System  *  Total  Pressure  Increase  Resulting  from  the  Bilateral 
Junction  as  Compared  with  the  Predicted  Total  Pressure  for  a  Two  Single-lateral  Setup 


Percent  Total  Pressure  Increase 

Velocitv  Ratio 

Averace 

Maximum 

0.3  -  3.7 

2.0 

3.5 

0.5  -  2.0 

1.7 

2.5 

0.75  -  1.25 

1.5 

2.3 

*  Pressure  measured  at  the  VP  measuring  station  in  duct  "m". 


When  interpreting  all  of  the  results  presented  above  it  is  important  to  remember  that  the 
predicted  power  losses  for  the  two  single-lateral  junctions  were  estimated  under  "ideal" 
conditions  and  that  any  deviation  from  these  conditions  would  increase  the  pressure 
requirement  for  the  double  single-lateral  setup.  In  addition,  the  single-lateral  junction  model 
by  Guffey  and  Fraser  was  fit  to  a  more  practical  range  of  velocity  ratios  and  then  applied  to 
very  extreme  velocity  ratios  in  this  research.  Therefore,  the  junction  power  loss  equation 
derived  for  the  dual  single-lateral  setup  may  underestimate  the  true  power  losses  at  these 
extreme  velocity  ratios.  This  probably  explains  the  divergence  between  the  measured 
junction  p>ower  losses  for  the  bilateral  junction  and  the  predicted  losses  for  the  dual  single¬ 
lateral  set-up  seen  at  extreme  velocity  ratios  in  Figure  47. 


CHAPTER  IV 
CONCLUSIONS 


The  linear  kinetic  power  model  provided  a  good  fit  to  the  empirical  data  collected  for  a 
bilateral  junction,  and  is  likely  to  produce  a  similarly  good  fit  for  any  other  bilateral 
junction. 

As  was  anticipated  from  momentum  considerations,  power  losses  were  least  when  the 
velocities  of  the  airstreams  in  the  converging  ducts  were  equal  and  were  greatest  when  the 
lateral  velocities  exceeded  the  velocity  of  the  center  duct.  Power  losses  increased  rapidly 
when  the  lateral  duct  velocities  exceeded  the  center  duct  velocity.  Powers  losses  increased 
to  a  lesser  degree  when  the  center  duct  velocity  exceeded  that  of  the  laterals. 

The  dynamic  power  losses  measured  in  the  bilateral  junction  were  14-92%  (average  of 
45%)  greater  than  the  power  losses  predicted  for  two  single-lateral  junctions  in  an 
equivalent  setup.  For  most  of  the  velocity  ratios  tested,  the  incremental  increase  in  power 
lost  for  the  bilateral  junction  was  less  than  the  power  loss  that  would  result  from  placing 
one  90°  elbow  (5  section,  R=2D)  or  an  additional  10  feet  of  duct  in  the  main.  For  some 
layouts,  avoiding  the  use  of  bilateral  junctions  may  require  additional  elbows  and/or 
additional  length  in  the  branch  ducts  as  shown  in  Figure  8.  In  these  cases,  a  bilateral 
junction  is  likely  to  be  more  efficient  than  two  single-lateral  junctions  in  series;  especially  if 
the  two  single  lateral  junctions  located  are  close  together. 

The  ASHRAE  table  of  coefficients  is  incorrect  and  should  not  be  used  to  design 
systems  containing  bilateral  junctions.  Designers  should  refer  to  Idelchik  for 
coefficients  for  bilateral  junction.  Caution  should  also  be  used  when  using  Idelchik's 
values  since  errors  were  also  found  in  this  table. 

The  linear  kinetic  power  model  was  effective  in  predicting  the  losses  in  the  bilateral 
junction  studied.  The  fit  to  data  was  acceptable  for  the  extreme  range  of  velocity  ratios 
tested  (0-3.7)  and  very  good  for  practical  velocity  ratios.  TTiC  linear  kinetic  power  model 
should  be  applied  to  empirical  data  from  a  diverse  set  of  bilateral  junctions  so  that  loss 
coefficients  can  be  modeled  as  functions  of  junction  geometry  as  has  been  done  for  single¬ 
lateral  junctions  Research  should  be  limited  to  a  velocity  ratio  range  of  0.5-2.0  since 
this  range  of  velocity  ratios  should  cover  most  design  situations. 
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For  many  situations,  a  design  incorporating  bilateral  junctions  may  have  a  lower  power 
requirement  than  one  which  msc;  two  single-latera'  junctions  in  place  of  each  bilateral 
junction.  Even  in  cases  where  the  single-lateral  setup  is  more  efficient,  the  advantages  of 
the  bilateral  junction  may  outweigh  the  small  potential  energy  savings  of  the  single-lateral 
junctions. 
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5-50  INDUSTHiAL  VENTILATION 


PREFERRED 


PREFERRED 


PREFERRED  ACCEPTABLE  AVOID 


BRANCH  ENTRY 

Branches  should  enter  at  gradual  expansions  and  at  on  angle 
of  DO° or  less(preferred)  to45°  if  necessary.  Expansion  should 
he  I5°maximum.  See  Fig.  6-12  for  Loss  Factor. 


PREFERRED  A  VOID 

proper  duct  SIZE 

Sire  the  duct  to  maintain  the  selected  or  higher 
transport  velocity. 


Appendix  B 

ASHRAE  Coefficients  for  Bilateral  Junctions 


Appendix  C 

Idelchik's  Coefficients  for  Bilateral  Junctions 


Four-way  wye  piece  of  the  type  f.j  =  F.j  =  =  F^. 

a  =  15°  [11.  12! 


Diagram 

7-25 


1.  Merging  of  streams  i  outlet  four-'.',  ay  piece) 
Side  branch 

^  =  1  - ' 


\QcF:,)  \QcJ 

[Qc  Q:s-<i  -^:s.-gi5ii’  ,  T  ( 

►Merging  '  ^  -  F,  ~  Q,„Q  \Qc  j  F..s^  V 

Dividing 

see  the  curves  ^  =  fyQ^i,  Qc>  Q'.pQvs^  at  different  F.^F^. 

For  the  other  side  branch  subscripts  1  and  2  change  places. 

Main  passage 

^  _  -^Pii  _  j  / 2st\  "  _  / ^si\  '  _ 1  ~  ^sr.  ^c 

uwj,2~  ^\Qc)  \Qcj  (0.75  -  0.25est-2cl‘ 

\Qc/  F-.F-'^  -Q:s,Q:s)-\Qsi  / 
see  the  curves  {•(.  jf  =  fiQst  Qc'  Q-.s'Q.s^  (*'  different  F.pF^.. 

2.  Division  of  now  lintake  four-way  piece),  j  and  are  determined  tentatively  similar  to  diverging 
lAyes  from  Diagrams  7-15  and  7-17  (curve  I) 


1 .0  0.5  ana  2.0 


Values  of  f. 


F.,F, 

=  0.2 

0,37 

0.46 

1.48 

0.29 

0.43 

1.23 

0.32 

-0.31 

-1.13 

F.^'F, 

=  0.4 

0.50 

-0.05 

0,34 

0.39 

0.06 

0.31 

0.27 

-0.10 

-0.65 

3 


101(0 


(.ContinueJ^  6 


Four-way  wye  piece  of  the  type 
a  =  15’  [11.  12) 

—  Pjs~  Ps- Pst  ~  Pc’ 

Diagram 

1-25 

Values  of  j, 

n  . 

Qst^Qc 

Q,,  0  Q.l  0.2 

0.3  0.4  Oi 

0.6 

0.7 

0.8  0.9  1.0 

FJF,  =  0.2 


0.5-2.0 

-4.37 

-2.93 

-2.04 

-1.44 

-1.08 

-0.58 

-0.22 

0.03 

0.16 

0.14 

0 

1.0 

-3.84 

-2.93 

-2.13 

-1.44 

-0.89 

-0.45 

-0.13 

0.08 

0.17 

0.14 

0 

PJPc  = 

0.4 

0.5 -2.0 

-1.70 

-1.19 

-0.76 

-0.40 

-0.12 

0.08 

0.21 

0.27 

0.25 

0.16 

0 

1.0 

-1.42 

-0.96 

-0.58 

-0.26 

-0.02 

0.15 

0.26 

0.29 

0.26 

0.16 

0 

PJPc  = 

0.6 

0.5-2.0 

-0.81 

-0.47 

-0.19 

0.04 

0.20 

0.30 

0.36 

0.35 

0.29 

0.17 

0 

1.0 

-0.16 

-0.31 

-0.05 

0.13 

0.27 

0.35 

0.39 

0.37 

0.29 

0.17 

0 

PjPc  = 

1.0 

0.5-2.0 

-0.35 

-0.11 

0.10 

0.26 

0.36 

0.42 

0.43 

0.39 

0.31 

0.18 

0 

1.0 

-0.21 

0.02 

0.19 

0.33 

0.41 

0.45 

0.45 

0.41 

0.31 

0.18 

0 

Diagram 

7-26 


Four-way  wye  piece  of  the  type  F^j  —  Pis  ~  Ps'<  Pst  ~  Pc< 
a=  30’  (11,  121 


see  the  curves  f,c.i  =  flQis/Qc-  Qsi/Q,s)  at  different /■,,//=■£.. 
For  the  other  side  branch  subscripts  1  and  2  change  places. 
Main  passage 


r  ^Pst  ■i./'gstV  fQstV  l^gsi/gc 

pw»/2  \QcJ  \Qc)  (0.75  +  0.25C5,/(?c)’ 

^  1  ^te.s/g.x)7gc  .y 

■  \Qc)  f.r  (1  +g.Vg.r)’  ^gst  J 

see  the  curves  f^.st  =AQst/Qc<  QjjQts)  a*  different  P,sJP? 

2)  Division  of  flow  (intake  four-way  piece);  j  and  jj  are  determined  tentatively  as  for  diverging  wyes 
from  Diagrams  7-15  and  7-17  (curve  1). 


374 


Four-way  wye  piece  of  the  type  F 
a  =  45°  Ill,  121 


yV2i.  F2x 


1)  Merging  of  streams  (outflow  four-way  piece) 
Side  t  ranch 


a  +<2,w'C.r)l 


_  I  42  ik  i±(QJQ^  Qc__  A  ’ 

■  \Qc)  (1  +  e,r/G.r)’  \Q^  ) 


see  the  curves  fc.st  =  RQulQo  Gia/Gir)  at  different  F^JF^. 

^c.s>  fc.st 

— — — ■  _ I - 1 — - 1  2)  Division  of  flow  (diverging  four-way  piece);  ^  and 

\p^efs  ~  Qj  I  fe.jt  are  determined  tentatively  as  for  wyes  from  Diagrams 

I  /  I  I  7-15  and  7-17  (curve  1). 


mm 


\^y/(  ho  aSinti  2.0  Values  of  . 


OA  \  0.6  \  Q.S 


^  C 


m 


z 


^-ij) 

'"c 


2.0  - 

0.5  - 

1.0  - 

2.0  - 


m 

SBSSZm 


PJFc 

=  0.2 

■■ 

■ 

-0.36  0J9 

1.77 

3.20 

4.88 

-0.24  0.63 

1.70 

2.64 

3.73 

-0.19  0.21 

0.04 

- 

- 

FjFc 

=  0.4 

0^8 

0.92 

1.31 

0.36  0.17 

0.55 

0.72 

0.78 

0.18  0.16 

-0.06 

- 

- 

=  0.6 

OJO  -0.07 

0.31 

0.60 

0.82 

0.37  0.12 

0.55 

0.60 

0.52 

0.18  0J6 

0.16 

- 

- 

l^JPc 

=  1.0 

0^1  -0.09 

0.25 

OJO 

0.65 

0.37  0.13 

0.46 

0.61 

0J4 

0.15  0.38 

0.42 

- 

3 
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Diagram 

7-27 


Values  of 


QjQc 


Q:s 


g,r 

0 

0.1 

0.2 

0.3  0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1. 

0. 5-2.0 

-2.92 

-1.87 

-1.29 

FjFc  =  0.2 
-0.80  -0.56 

-0.23 

-0.01 

0.16 

0.22 

0.15 

0 

1.0 

-2.54 

-1.87 

-1.30 

-0.80  -0.42 

-0.12 

0.08 

0.20 

0.22 

0.15 

0 

0.5-2.0 

-0.98 

-0.61 

-0.30 

FJF,  =  0A 
-0.05  0.14 

0.26 

0.33 

0.34 

0.28 

0.17 

0 

1.0 

-0.77 

-0.44 

-0.16 

0.05  0.21 

0.31 

0.36 

0.35 

0.29 

0.17 

0 

0.5-2.0 

-0.32 

0.08 

0.11 

F.VFc  =  0.6 
0.27  0.37 

0.43 

0.44 

0.40 

0.31 

0.18 

0 

1.0 

-0.18 

-0.04 

0.21 

0.34  0.42 

0.46 

0.46 

0.41 

0.31 

0.18 

0 

0.5-2.0 

0.11 

0.36 

0.46 

FJF,  =  1.0 
0.53  0J7 

0.56 

0J2 

0.44 

0.33 

0.18 

0 

1.0 

0.29 

0.42 

0.51 

0.57  0.59 

0.58 

0J4 

0.45 

0.33 

0.18 

0 

Four-way  wye  piece  of  the  type  F,  ^  —  F,  j  =  Ff,  F51  —  Fj ;  Diagram 

a  =  60*  (11.  12)  ■  7-28 


1)  Merging  of  streams  (converging  four-way  piece). 
Side  branch 


-8 


(Qj\'  [Qc/0.r-(l  ^g,s/g.r)l’ 
\Qc/ 


4-(i  *QjQ,s)QjQc 


[s. 


1  -I- 


see  the  curves  =  fiQjQc<  QijQy^  at  different  F^JF^. 

For  the  other  side  branch  subscripts  1  and  2  change  places. 

Main  passage 

j.  ^Pst  _  .  .  / gsA  ’  _  / gjA  *  1  g$t/gc _ / gsA  ’  ^ 

\Qc)  \Qc)  (OJS +0.2SQJQc)  \QcJ  F^ 


1  +  (g.^g.r)’ 

[1  +(g,r/g,s)|’ 

Vg.t 

see  the  curve  =  flQulQc,  at  different  F.j/Fg. 

2)  Division  of  flow  (diverging  four-way  piece);  f<..jand  are  determined  tentatively  as  for  diverging 
wyes  from  Diagram  7-17  (curve  1). 


I 


(ContinuedyQ  Q 


Four-way  wye  piece  of  the  type  Diagram 

a  =  60“  [11,12)  7-28 


?C.J '  fc-St 


Values  of  ^,..1 

QjQc 

QiS  - 


Q,s 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.5 

-1.0 

-0.31 

0.59 

=  0.2 

2.00 

3.62 

5.54 

7.72 

1.0 

-1.0 

-0.20 

0.80 

2.07 

3.30 

4.77 

- 

2.0 

-1.0 

-0.09 

0.62 

0.97 

- 

- 

- 

0.5 

-1.0 

-0.47 

-0.06 

=  0.4 

0.60 

1.12 

1.63 

2.10 

1.0 

-1.0 

-0.34 

0.25 

0.73 

1.10 

1.31 

- 

2.0 

-1.0 

-0.15 

0.27 

0.41 

- 

- 

- 

0.5 

-1.0 

-0.50 

0.04 

=  0.6 

0.38 

0.74 

1.03 

1.23 

1.0 

-1.0 

-0.36 

0.18 

0.67 

0.82 

0.87 

- 

2.0 

-1.0 

-0.15 

0.40 

0.47 

- 

- 

- 

0.5 

-1.0 

-OJO 

FJF, 

-0.07 

=  1.0 

0.30 

0.58 

0.79 

0.88 

1.0 

-1.0 

-0.36 

0.16 

0J3 

0.74 

0.75 

- 

2.0 

-1.0 

-0.13 

0.46 

0.61 

- 

- 

- 

Values  of 


Qxx/Qc 


Wax 

G.r 

0 

0.1 

0.2 

0.3  0.4 

0.5 

0.6 

0.7 

0  j 

0.9 

1.0 

OJ-2.0 

-1.77 

-1.02 

-0.64 

FJF,  =  0.2 
-0.30  -0.15 

0.06 

0.20 

0.26 

0.26 

0.16 

0 

1.0 

-IJO 

-1.03 

-0.64 

-0.30  -0.05 

0.13 

0.24 

0.29 

0.26 

0.16 

0 

0.5-2.0 

-0.40 

-0.14 

0.07 

FJF,  =  0A 
0.24  0.35 

0.41 

0.42 

0.39 

0.30 

0.18 

0 

1.0 

-0.25 

-0.02 

0.16 

0.31  0.40 

0.44 

0.45 

0.40 

0.31 

0.18 

0 

0.5-2.0 

0.06 

0.23 

0.36 

FjFc  =  0.6 
0.46  0 J I 

0.52 

050 

0.43 

0.32 

0.18 

0 

1.0 

0.16 

0.32 

0.43 

0.51  0.55 

0.55 

051 

0.44 

0.33 

0.18 

0 

OJ-2.0 

0.44 

0J4 

0.60 

FjFc 

0.65  0.65 

0.62 

056 

0.47 

0.34 

0.18 

0 

1.0 

OJO 

0.59 

0.64 

0.67  0.67 

0.63 

057 

0.47 

0.34 

0.18 

0 

378 


Cross  ot  the  type  F 


F.;a  =  90‘' 


Diagram 


*  Merging 
-  Dividing 


W2^F2s  i 
Main  passage 

fc.i'  fc.st 


iQa'Qis-n  ^-g=/g,s)r 

\Qc  J  (I  Qis/'Q is^O is'^c 

see  the  curves  fic.^  = /(Qis/Qc.  2:r>'Qis)  different  F^jF^. 
For  the  other  side  branch  subscripts  1  and  2  change  places. 


.  _  ^Pst 

■  pw*/2 


I  i  !  I  ! 

^':c£S,f£^3riciF.0\ 

T'^J  !  I 


•^Pst  _  J  I  ^st\  '  ( ^StN  ’ _ I  QixlQc _ 

\Qc)  V  2cy  (0.75  ^  0.25(2st/2c)‘ 

see  the  curves  =  iXQ^JQc,  QiJQis)  a'  different  F^  JF^ 
For  standard  crosses  made  of  malleable  iron  at  Q^lQc  ^  0.7 

2)  Flow  division  (diverging  cross);  f ^  and  jj  are  deter¬ 
mined  tentatively  as  for  diverging  wyes  from  Diagrams  7-15 
and  7-17  (curve  1) 

Values  of 


Bnrij 

UWil 


O.Sii03nf\2.0 


at 

a 

'6 

1 

iEiBni 

ninn 


FjFc 

=  0.2 

-0.10 

1.09 

2.72 

4.77 

7.25 

-0.05 

1.35 

3.12 

5.05 

7.40 

-0.31 

1.77 

3.37 

- 

- 

FJF, 

=  0.4 

-0.29 

0.34 

1.03 

1.77 

2.56 

■0.14 

0.60 

1.33 

2.05 

2.71 

0.12 

1.02 

1.68 

- 

- 

FJPc 

=  0.6 

•0.32 

0.20 

0.72 

1.22 

1.70 

■0.18 

0.46 

1.02 

IJO 

1.85 

0.09 

0.88 

1.37 

- 

- 

=  0.8 

0.33 

0.13 

0.61 

1.02 

1.38 

0.18 

0.41 

0.91 

1.30 

1.54 

0.08 

0.83 

1.26 

- 

- 

FJF^ 

=  1.0 

0.34 

0.13 

0.56 

0.93 

1.25 

0.19 

0,39 

0.86 

1.21 

1.40 

0.07 

0.81 

1.21 

- 

- 

0  01  0-7  O  f  0.4  0.i  0.6  0.7  0.8  0.9  1.0 


1.20  1.19  1.17  1.12  1.03  0.96  0.85  0.72  0.56  0.39  0  20 
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Appendix  D 

Corrected  Coefficients  for  ASHRAE  Table  5-34  based  on  Formulas 

Provided  by  Idelchik 
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Coefficients  for  Laterals 


Qa/Qm 


QbZQa 

a 

ILl 

iLI 

1L2  0.4 

0.5 

Aa/Am=0.2 

0.5 

7 

-0.34 

0.59 

1.77  3.20 

4.88 

1 

7 

-0.24 

0.63 

1.60  2.63 

3.70 

2 

7 

-0.16 

0.20 

0.03  -0.79 

-2.43 

Aa/Am=0.4 

0.5 

7 

-0.48 

0.01 

0.48  0.91 

1.30 

1 

7 

-0.36 

0.17 

0.55  0.76 

0.79 

2 

7 

-0.17 

0.16 

-0.06  -0.95 

-2.68 

Aa/Am=0.6 

0.5 

7 

-0.50 

-0.07 

0.30  0.59 

0.80 

1 

7 

-0.37 

0.12 

0.45  0.59 

0.51 

2 

7 

-0.15 

0.26 

0.16  -0.56 

-2.06 

Aa/Am=0.1 

0.5 

7 

-0.51 

-0.09 

0.25  0.50 

0.65 

1 

7 

-0.37 

0.13 

0.46  0.61 

0.54 

2 

7 

-0.12 

0.38 

0.43  -0.09 

-1.33 

Notes: 

a 

=  lateral 

b 

=  lateral 

m  =  downstream  main 
?=  Equation  has  a  zero  in  the  denominator, 
source  of  table  values  is  uknown. 


0.6 


6.78 

4.77 

-5.11 


1.63 

0.58 

-5.47 


0.91 

0.18 

-4.59 


0.70 

0.22 

-3.52 
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Notes:  a  =  lateral 

b  =  lateral 

m  =  downstream  main 
s  =  collinear  duct 

?=  Equation  has  a  zero  in  the  denominator, 
source  of  table  values  is  uknown. 


Appendix  E 

Calibration  Data  Worksheet. 


Appendix  F 

Equipment  Calibration  Curves 


Calibration  Curve  Summary 
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Inclined  Manometer  #  1 

Pressure  =  0.998487  (Inc#l  reading)  -  0.002390 

Inclined  Manometer  #  2 

Pressure  =  0.995042  (Inc#2  reading)  -  0.003129 

Merical 

Pressure  =  0.997269  (Merical  reading)  +  0.005857 

Transducer  #  5 

Pressure  =  1.89691  (Transducer  #5  Volts) 

+  0.009559  (Transducer  #5  Volts)^-  6.70500 


Transducer  #  3 


Pressure  =  7.81296  (Transducer  #3  Volts) 

-  0.033744(Transducer  #3  Volts)^  -  0.769084 


Transducer  #  4 


Pressure  =  1.98680  (Transducer  #4  Volts)  -  7.00087 


Transducer  #  8 


Pressure  =  1.79956  (Transducer  #8  Volts) 

+  0.020732(Transducer  #8  Volts)^  -  6.63 1 19 


Regression  Ression  Results 
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Inclined  Manometer  #1 


Dependent  variable  is:  Corrected  Hook  Gage 

446  total  cases  of  which  141  are  missing 
RD  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0041  with  305  -  2  =  303  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  364.409 

Residual  0.005183 


df  Mean  Square  F-ratio 

1  364  21303754 

303  0.000017 


Variable 

Constant 

Incl 


Coefficient 

s.e.  of  Coeff 

t-ratio 

-0.002390 

0.0004 

-6.25 

0.998487 

0.0002 

4616 

Pressure  =  0.998487  (Inc#l  reading)  -  0.002390 


Inclined  Manometer  #2 

Dependent  variable  is:  Corrected  Hook  Gage 

446  total  cases  of  which  141  are  missing 
Rn  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0040  with  305  -  2  =  303  degrees  of  freedom 


Source 

Sum  of  Squares  df 

Mean  Square 

F-ratio 

Regression 

364.410 

1 

364 

22921968 

Residual 

0.004817 

303 

0.000016 

Variable 

Coefficient 

s.e.  of  Coeff 

t-ratio 

Constant 

-0.003129 

0.0004 

-8.48 

Inc2 

0.995042 

0.0002 

4788 

Pressure  = 

0.995042  (Inc#2  reading) 

-  0.003129 

Merical 


1 


Dependent  variable  is:  Corrected  Hook  Gage 

446  total  cases  of  which  309  are  missing 
Rn  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0069  with  137  -  2  =  135  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  1590.69 

Residual  0.006429 


df 

Mean  Square 

F-ratio 

1 

1591 

33401107 

135 

0.000048 

Variable 

Constant 

Merical 


Coefficient 

s.e.  of  Coeff 

t-ratio 

0.005857 

0.0009 

6.34 

0.997269 

0.0002 

5779 

Pressure  =  0.997269  (Merical  reading)  +  0.005857 


Transducer  #5 

Dependent  variable  is:  Corrected  Hook  Gage 

cases  selected  according  To  HKcorr>0 
446  total  cases  of  which  130  are  missing 
Ra  =  100.0%  R2(adjusted)  =  100.0% 
s=  0.0077  with  316-3  =  313  degrees  of  freedom 


Source 

Sum  of  Squares  df 

Mean  Square 

F-ratio 

Regression 

619.751 

2 

310 

5249846 

Residual 

0.018475 

13 

0.000059 

Variable 

Coefficient 

s.e.  of  Coeff 

t-ratio 

Constant 

-6.70500 

0.0156 

-430 

Ch5 

1.89691 

0.0068 

280 

Ch5SQ 

0.009559 

0.0007 

13.3 

Pressure  =  1.89691  (Transducer  #5  Volts) 

+  0.009559  (Transducer  #5  Volts)^  -  6.70500 


Transducer  #  3 
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Dependent  variable  is:  Corrected  Hook  Gage 

cases  selected  according  To  HKcorr>0 
446  total  cases  of  which  194  are  missing 
Ro  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0055  with  252  -  3  =  249  degrees  of  freedom 


Source 

Sum  of  Squares  df 

Mean  Square  F-ratio 

Regression 

316.568 

2  158 

5262728 

Residual 

0.007489 

249 

0.000030 

Variable 

Coefficient 

s.e.  of  Coeff 

t-ratio 

Constant 

-0.769084 

0.0017 

-451 

Ch3 

7.81296 

0.0118 

664 

Ch3SQ 

-0.033744 

0.0175 

-1.93 

Pressure  =  7.81296  (Transducer  #3  Volts) 

-  0.033744(Transducer  #3  Volts)2  -  0.769084 


Transducer  #4 


Dependent  variable  is:  Corrected  Hook  Gage 

cases  selected  according  To  HKcorr>0 
446  total  cases  of  which  151  are  missing 
Ra  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0084  with  295  -  2  =  293  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  344.584  1 

Residual  0.020569 


df  Mean  Square  F-ratio 
345  4908527 

293  0.000070 


Variable 

Constant 

Ch4 


Coefficient  s.e.  of  Coeff  t-ratio 

-7.00087  0.0038  -1825 

1.98680  0.0009  2216 


Pressure  =  1.98680  (Transducer  #4  Volts)  -  7.00087 


Transducer  #8 

Dependent  variable  is:  Corrected  Hook  Gage 

cases  selected  according  To  HKcoit>0 
446  total  cases  of  which  156  are  missing 
Ro  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0102  with  290  -  3  =  287  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  337.586  2 

Residual  0.029574 


df  Mean  Square  F-ratio 
169  1638047 

287  0.000103 


Variable 

Constant 

Ch8 

Ch8sq 


Coefficient 

s.e.  of  Coeff 

t-ratio 

-6.63119 

0.0388 

-171 

1.79956 

0.0177 

102 

0.020732 

0.0020 

10.4 

Pressure  =  1.79956  (Transducer  #8  Volts) 

+  0.0207 32(7 ransducer  #8  VoJts)2-  6.63119 


Appendix  G 
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Duct  -  20  and  24  gage  spiral  wound  duct,  connectors,  expansions,  and  elbows  purchased 
from  United  McGill  Corporation  .Stockton  CA. 

Junction  -  Bilateral  junction  with  5"  diameter  side  laterals  entering  at  a  45°  angle  from 

centerline;  7"  diameter  upstream  coaxial  duct;  and  10"  diameter  downstream  main. 
Expanded  from  7"  to  10"  diameter  over  7".  Junction  was  purchased  from  United 
McGill  Corporation,  Stockton  CA. 

Fan  -  Buffalo  Forge  MW40  Centrifugal  Fan,  Buffalo  NY. 

Motor  -  Lincoln  Motors,  AC,  20  HP 

Transmission  -  Infinitely  Variable  Mechanical  Transmission  Model  7-10A(E)  0:1-1  CCW 
by  Speed  Control  Industries  Inc,  Richland  WA. 

Electronic  Tachometer  -  DITAK  5  Model  DT-5  by  Red  Lion  Controls 

Duct  Tape 

Pitot  Tubes  -  Dwyer  hemispherical  probe  Pitot  tubes,  3/16" 

Pitot  Tube  Holders  -  Guffey 

Tygon®Tubing  -  American  Scientific  Products  T6000-4,  3/16  ID  x  1/16  wall  thickness 

High  density  polypropylene  quick  disconnect  fittings 

Inclined  Manometer  -  Meriam  Instrument,  Model  40HE35WM 

Hook  Gage  -  Dwyer  Instruments,  Inc.,  Series  1425 

Indicating  Fluids: 

Meriam  1000  Green  Concentrate,  No.  922WA 

Dwyer  Fluorescein  Green  Concentrate,  1.000  specific  gravity 


Hand  Pump  -  Meriam  Instruments,  Model  B34348 
Pressure  Transducers: 

Omega  Engineering,  Inc.,  PX160  and  PX  154  Series 
Meri-cal®  Model  DP200I,  Meriam  Instruments 
Transducer  Power  Supply  -  Omega®  PST-4I30  regulated  power  supply. 

Barometer  (Mercurial)  -  Princo  Instruments,  Inc.,  No.  453 
Thermometers: 

Glass  Thermometer  -  Ever  Ready  Thermometer  Company,  Instrument  No.  1446, 
ASTM  Precision  Specifications 

Digital  -  Fluke  52  K/J  Thermometer 

Psychrometer  -  Cole-Parmer  Instrument  Co.,  Psychro-Dyne,  No.  3312-40 

Measuring  Tapes: 

Ace  Hardware,  No.  24988,  25  ft. 

Empire  Model  6550,  50  ft. 

Stanley  “Powerlock“,  No.  33-203,  3  m/10  ft. 

Caliper  (made  in  China  -  manufacturer  unknown) 

Level  -  construction-type 
Macintosh  SE  Computer 

Thumb  Wheel  -  Cherry  Programmable  Controller  Peripheral,  TH  Series 
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Computer  Programs: 

Analog  Connection  Workbench™  2.0  by  Strawberry  Tree  Computers  Inc. 
Microsoft®  Excel,  Microsoft  Corporation 
Data  Desk®3.0,  Odesta  Corportation 


Appendix  H 

Friction  Data  Collection  Procedure 
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Friction  Data  Collection  Procedure 

1 .  Calibrate  equipment. 

2.  Clean  duct. 

3.  Assemble  and  seal  duct. 

4.  Drill  holes  in  duct  for  differential  static  pressures. 

5.  Install  Pitot  tubes  and  Pitot  tube  holding  devices. 

6.  Conduct  leak  test  on  duct  and  Tygon®  tubing. 

7.  Check  alignment  of  duct. 

8.  Start  "Workbench"  program  and  insure  transducer  is  on  for  at  least  30  minutes 

(for  warm-up)  before  collecting  data. 

9.  Turn  on  Fluke  and  allow  to  warm  up  for  at  least  15  minutes. 

10.  Check  Fluke  against  mercury  thermometer  and  adjust  if  necessary. 

11.  Zero  barometer  according  to  manufacturers  instructions. 

12.  Turn  fan  on  and  set  rpm  to  achieve  the  desired  duct  velocity. 

13.  Take  the  following  measurements  for  the  measurement  station. 

a.  Record  "Run"  number  on  computer. 

b.  Record  fan  rpm  on  data  collection  worksheet  (check  reading  periodically 

and  note  range). 

c.  Measure  airstream  temperature  just  upstream  of  the  measuring  station  and 

record  on  data  collection  worksheet. 

d.  Conduct  first  10  point  Pitot  traverse  and  collect  centerline  VP  on 

transducer  (BCD  1-10  and  31). 

e.  Conduct  second  10  point  Pitot  traverse  and  collect  centerline  VP  on 

transducer  (BCD  11-10  and  32). 

f.  Conduct  third  10  point  Pitot  traverse  and  collect  centerline  VP  on 

transducer  (BCD  21-30  and  33). 

g.  Record  centerline  static  pressure  of  third  traverse  with  transducer 

(BCD  34). 

h.  Check  level  and  zero  of  inclined  manometer. 
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i.  Measure  centerline  static  pressure  of  third  traverse  with  inclined 

manometer  and  record  on  data  collection  worksheet. 

j.  Measure  centerline  velocity  pressure  of  third  traverse  with  inclined 

manometer  and  record  on  data  collection  worksheet. 

k.  Measure  wet  and  dry  bulb  temperatures  and  the  barometric  pressure 

and  record  on  the  data  collection  worksheet. 

14.  Record  the  transducer  voltages  for  all  static  pressure  differentials  to  a  computer 

file. 

15.  Recheck  the  level  and  zero  on  the  inclined  manometer. 

16.  Measure  all  static  pressure  differentials  with  an  inclined  manometer  and  record 

the  data  on  the  data  collection  worksheet. 

17.  Adjust  the  fan  rpm  to  achieve  a  new  duct  velocity  and  return  to  step  13. 

18.  Repeat  until  adequate  data  points  are  taken  within  the  velocity  range  of  interest. 

19.  Change  setup  to  a  new  duct  diameter  and  repeat  steps  2  through  18. 


Appendix  I 

Friction  Data  Collection  Worksheet 


Friction  Data  Collection  Worksheet  Pag0 


da/  SP  6-6 
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Junction  Data  Collection  Procedure 
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Junction  Data  Collection  Procedure 

1 .  Calibrate  equipment. 

2.  Clean  duct. 

3.  Complete  friction  estimation  for  ducts. 

4.  Grind  inside  of  junction  to  remove  excess  metal. 

5.  Connect  ducts  to  junction  and  install  velocity  pressure  measuring  stations. 

6.  Drill  holes  to  measure  hood  static  pressures  and  airstream  temperatures. 

7.  Check  level  of  duct,  straighten  with  stringline  and  seal  with  duct  tape.  Seal  measuring 

station  connections  with  duct  tape  and  silicone  sealant. 

8.  Install  Pitot  tubes  and  Pitot  tube  holders  and  zero. 

9.  Conduct  leak  test  on  system  and  Tygon®  tubing. 

10.  Start  "Workbench"  program  and  insure  transducer  is  on  for  at  least  30  minutes  (for 

warm-up)  before  collecting  data. 

1 1.  Turn  on  Huke  and  allow  to  warm  up  for  at  least  15  minutes. 

12.  Check  Fluke  against  mercury  thermometer  and  adjust  if  necessary. 

13.  Zero  barometer  according  to  manufacturers  instructions. 

14.  Set  desired  branch  velocity  ratios  by  adjusting  the  blast  gates  in  each  branch  using 

impedance  method  Adjust  final  velocities  by  using  the  variable  drive  transmission. 

15.  Check  the  level  and  zero  of  the  inclined  manometer  and  adjust  as  necessary. 

16.  Move  control  box  to  VP  measuring  station  in  the  main. 

17.  Connect  transducer  and  inclined  manometer  to  the  hand  pump  with  Tygon®  tubing  and 

put  under  pressure. 

18.  Log  calibration  pressure  on  the  computer.  Record  inclined  manometer  reading  on  data 
collection  worksheet. 

19.  Take  measurements  for  duct  M. 

a.  Record  "Run"  number  on  computer. 

b.  Record  fan  rpm  on  data  collection  worksheet  (check  reading  periodically  and 

note  range). 

c.  Measure  airstream  temperature  just  upstream  of  the  measuring  station  and 

record  on  data  collection  worksheet. 
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d.  Conduct  first  10  point  Pitot  traverse  and  collect  centerline  VP  on  transducer 

(BCD  1-10  and  31). 

e.  Conduct  second  10  point  Pitot  traverse  and  collect  centerline  VP  on  transducer 

(BCD  11-10  and  32). 

f.  Conduct  third  10  point  Pitot  traverse  and  collect  centerline  VP  on  transducer 

(BCD  21-30  and  33). 

g.  Record  centerline  static  pressure  of  third  traverse  with  transducer  (BCD  34). 

h.  Measure  static  pressure  differential  between  the  centerline  of  M  and  the 

centerline  of  S  with  the  transducer  (BCD  35). 

i.  Check  level  and  zero  of  inclined  manometer. 

j.  Measure  static  pressure  differential  between  the  centerline  of  M  and  the 

centerline  of  S  with  inclined  manometer  and  record  on  data  collection 
worksheet. 

k.  Measure  centerline  static  pressure  of  third  traverse  with  inclined  manometer 

and  record  on  data  collection  worksheet. 

l.  Measure  centerline  velocity  pressure  of  third  traverse  with  inclined 

manometer  and  record  on  data  collection  worksheet. 

m.  Measure  wet  and  dry  bulb  temperatures  and  the  barometric  pressure 

and  record  on  the  data  collection  worksheet. 

20  Move  the  control  box  to  the  junction. 

21  Check  the  level  and  zero  of  inclined  manometer. 

22.  Connect  transducer  and  inclined  manometer  to  the  hand  pump  with  Tygon®  Uibing 

and  put  under  pressure. 

23.  Log  calibration  pressure  with  the  control  box  at  the  junction.  Record  inclined 

manometer  reading  to  data  collection  worksheet. 

24.  Take  measurements  for  duct  S  (repeat  steps  19  a  through  19  m  above). 

25.  Take  measurements  for  duct  A  (repeat  steps  19  a  through  19  m  above  except 

measure  the  static  pressure  differential  between  A  and  S). 

26.  Take  measurements  for  duct  B  (repeat  steps  19  a  through  19  m  above  except 

measure  the  static  pressure  differential  l^tween  B  and  S). 

27.  Return  to  step  14  above  for  new  setup. 

28.  Return  to  step  10  at  the  start  of  a  new  day. 


Appendix  K 
Impedance  Method 
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Impedance  Method 

The  ratio  of  the  airflows  for  each  upstream  branch  was  set  using  the  impedance  method 
(Guffey,  unpublished  manuscript)  and  the  "Connect"  and  "Balance"  programs  developed 
by  Dr.  Guffey.  To  use  the  impedance  method  a  baseline  characterization  of  the  system  was 
first  accomplished.  The  baseline  was  accomplished  with  all  blast  gates  fully  open  and  the 
fan  at  1130  revolutions  per  minute  (rpm).  The  static  pressure  just  upstream  of  the  junction 
and  the  hood  static  pressure  (approximately  5  diameters  upstream  of  the  blastgate)  were 
measured  for  each  branch.  The  flow  in  each  branch  was  estimated  using  the  following 
equation: 


Average  Velocity  =  4005  V  0.85  VPcenterline 

The  physical  layout  of  system  components  were  entered  into  the  "Connect"  program. 
The  resulting  file  was  opened  by  the  "Balance"  program  and  the  baseline  system 
characteristics  were  entered.  Once  the  baseline  characteristics  were  entered,  the  target  flow 
for  each  branch  was  entered.  As  output,  the  program  gave  the  ratio  of  junction  static 
pressure  (SPJ)  to  hood  static  pressure  (SPH)  required  for  each  branch  to  achieve  the 
desired  flow.  Next,  the  blast  gate  in  each  branch  was  adjusted  until  the  desired  SPJ/SPH 
was  achieved.  The  centerline  velocity  for  each  branch  was  checked  to  ensure  the  desired 
velocity  ratios  were  attained.  If  required,  the  dampers  were  adjusted  slightly  to  fine-tune  the 
velocity  ratios.  The  blast  gates  were  fixed  in  position  once  the  desired  velocity  ratios  were 
achieved.  The  fan  rotation  rate  was  then  adjusted  to  attain  the  desired  velocities. 

The  Impedance  Method  of  balancing  the  system  worked  very  well  considering  the 
imprecision  in  the  method  of  estimating  the  baseline  air  flows. 
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Junction  Data  Collection  Worksheet 


Pressure  (BCD  -  78} 
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Summary  of  Empirical  Equations  For  All  Ducts  for  All  Test 

Sections 


4  inch  diameter 

/V  \  1-851  /  \ 

SPF4"  =  [TOO'}  fo*"  N.T.P.  in  4"  spiral,  fully  developed  flow 

s.e.  =  0.02457,  df  =  1 1 
SP12  =  1.0325  *  SP14 
SP13  =  1.0015  *  SP14 
SP23  =  0.984  *  SP14 
SP56  =  40.5"  =  1.85  *  SP14 


5  inch  diameter 

/V  \  1-85843  y  I 

SPF5"  =  0.3268  (lO^  }  N.T.P.  in  5"  spiral,  fully  developed  flow 

s.e.  =  0.0058,  df=  10 
SP12  =  1.0352  *  SP14 
SP13  =  0.97205  *  SP14 
SP23=  n/a 
SP34  =  1.083  *SP14 
SP56  =  1.6781  *  SP14 


6  inch  diameter 


SPF6"  = 

0.2854  j 

['.Yu*. 'I 
Uooo; 

1.8342 

I 

f  Li4  ' 
100' 

j,  for  N.T.P.  in  6"  spiral,  fully  developed  flow 

SPF6"  = 

0.2766  j 

uoooj 

1.7502 

1 

\ 

/  Li2  ■ 
1^100* 

s.e.  =  0.01 106" 

^  for  N.T.P.  in  6"  spiral,  fully  developed  flow 

SPF6"  = 

0.2468  ( 

UOOOy 

^  Li3  ^ 
100'} 

,  for  N.T.P.  in  6"  spiral,  fully  developed  flow 

SPF6"  = 

0.2774  \ 

uooo> 

,.s„( 

lOO'j 

,  for  N.T.P.  in  6"  spiral,  fully  developed  flow 

SPF6”  = 

0.3615  { 

uooo> 

'  L56  ^ 
lOO'j 

,  for  N.T.P.  in  6"  spiral,  fully  developed  flow 

SP12  =  0.96285  *  SP14 
SP13  =  0.99403  *  SP14 
SP23=  n/a 

SP34=  1.0603  *SP14 
SP56  =  1.375  *  SP14 
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SPF7"L14 

flow 


SPFIO"  = 


7  inch  diameter 


0.2156  N.T.P.  in  7"  spiral,  fully  developed 

SPL7"12/100  =0.224100*  (V/ 1000)^  1.847 3 

s.e.  =  0.01073  L14 

SP12  =  1.0433  *  SP14  -same  exponent 

SP13  =0.9975  *  SP14 

SP23  =  0.9502  *  SP14 

SP56  =  1.2767  *  SP14 


10  inch  diameter 

/V  \  1-889  /' J.,  \ 

(.1^)  (TO^  r  N-T.P.  in  10"  spiral,  fully  developed  flow 

s.e.  =  0.01048" 

SP12  =  1.0118  *  SP14  :  (mSP151/L12)/(uSP14/L14) 


SP23=  1.0578  *SP14  :(mSP154/L23)/(uSP14/L14) 


Regression  Analyses  for  Friction 


SPF  6" 

SPF  6'7100'  =  10'^-6.60445  *  V^l.85735  *  1007(447.6'712")  '  whole  range 
SPF  67100’  =  10^-6.61150  *  V^l.8600  *  100'/(447.6"/12”)  '  1000<V3<5000 
SPF  6"/100’  =  0.2551  *  (73/1000)^1.840 

fggfgssion  below  ***************** 

data  collected  9- 1 1  Jan90,  re-analyzed  30  March  90 

note:  must  do  9-10  WB  file  once  for  each  CF  file 
comment:  zero  shift  on  tSP,  especially  tSP12 

SP14/100'  =  0.2854  *  (V/1000)''l. 83417 
SP13/100'  =  0.246787*(V/1000)''1.8339 
SP12/100'=0.276589*(V/1000)''1.75016 
SPF34/100’  =  0.277487*(V/1000)''1.8127 
SPF56/100  =  0.361497*(V/1000)''l. 79686 

'  *****  show  that  transducer  predicts  VP  well  ******* 

Dependent  variable  is:  mVPcl 

24  total  cases  of  which  1  are  missing 
R2  =  100.0%  R2(adjustcd)  =  100.0% 
s=  0.0093  with  23-2  =  21  degrees  of  freedom 


Source 

Sum  of  Squares  df 

Mean  Square 

F-ratio 

Regression 

23.5497 

1 

23.55 

269696 

Residual 

0.001834 

21 

0.000087 

Variable 

Coefficient 

s.e.  of  Coeff 

t-ratio 

Constant 

0.003403 

0.0033 

1.04 

tVPclAvg 

0.998343 

0.0019 

519 

SPF  6”(cont) 


SP14  *************** 
Summary  statistics  for  PrErr  for  cmSP14 

cases  selected  according  To  selV  1500<V<60(X) 
NumNumeric  =17 
Mean  =  -0.00043 
Standard  Deviation  =  0.01355 


Dependent  variable  is:  logcmSP  1 4 

case*' selected  according  To  selV  1500<V<6000 
24  total  cases  of  which  7  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0057  with  17  -  2  =  15  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  1.33947 

Residual  0.000483 


df  Mean  Square  F-ratio 

1  1.3395  41604 

15  0.000032 


Variable 

Constant 

logVn 


Coefficient  s.e.  of  Coeff  t-ratio 
-1.01577  0.0053  -192 

1.83417  0.0090  204 


SP14  for  447.6"  =  10^-1.01577*(V/1000)^1.83417 
0.096434 


5p^"  12  ************** 

Dependent  variable  is:  logcmSP  12 

cases  selected  according  To  selV  2000<V<6000 
24  total  cases  of  which  9  are  missing 
R2  =  99.8%  R2(adjusted)  =  99.8% 
s  =  0.0083  with  15-2=13  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  0.555646 
Residual  0.000889 


df  Mean  Square  F-ratio 

1  0.555646  8126 

13  0.000068 


Variable 

Constant 

logVn 


Coefficient  s.e.  of  Coeff  t-ratio 

-1.64391  0.0120  -137 

1.75016  0.0194  90.1 


SPl  2=  10^-1. 6439 1*(V/1000)'' 1.7501 6 
SP 1 2=0.022703*(V/1000)^  1 .75016 
SP12/100'=10^-1.64391*100/(98.5/12)*(V/1000)^l. 75016 
SP12/100'=0.276589*(V/1000)''1.75016 


SPF  6"  (cent) 


§P|3  ***************** 

Summary  statistics  for  PrErr 

cases  selected  according  To  selV 
NumNumeric  =  22 
Mean  =  -0.00195 
Standard  Deviation  =  0.01055 


Dependent  variable  is:  logcmSPlS 

cases  selected  according  To  selV  2000<V<6000 
24  total  cases  of  which  2  are  missing 
r2  =  99.9%  R^fadjusted)  =  99.9% 
s  =  0.01 10  with  22  -  2  =  20  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  3.87226 
Residual  0.002438 


df  Mean  Square  F-ratio 

1  3.8723  31764 

20  0.000122 


Variable 

Constant 

logVn 


Coefficient  s.e.  of  Coeff  t-ratio 
-1.34419  0.0059  -228 

1.83392  0.0103  178 


SP13=0.045270*(V/1000)''l. 83392 

SP13/100'  =  0.045270*  100/(220. 125/12)*(V/1000)^1. 83392 

SP13/100'  =  0.24678/*(V/l  000)^1. 8339 


’***********  §p^"  34  ***>•<********** 

Summary  statistics  for  PrErr 

cases  selected  according  To  selV 
NumNumeric  =  16 
Mean  =  0.(XX)15 
Standard  Deviation  =  0.(X)690 


Dependent  variable  is:  logcmSP34 

cases  selected  according  To  selV  15{X)<V<6(XX) 
24  total  cases  of  which  8  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s=  0.(X)63  with  16-2  =  14  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  1.30795 

Residual  0.000558 


df  Mean  Square  F-ratio 

1  1.3079  32818 

14  0.000040 


Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constan  1.81270  0.0100  181 


SP6"=0.052607*(V/1 000)^  1 .8 1 27 

SPF67100'  =  0.052607*  100/(227.5/12)*(V/1000ri. 8 127 

SPF67100'  =  0.277487*(V/1000)''1.8127 


SPF  6"(cont) 


'*************  §p55  ************ 

Dependent  variable  is:  logcmSP56 

cases  selected  according  To  selV 
24  total  cases  of  which  9  are  missing 
r2  =  99.7%  R2(adjusted)  =  99.7% 
s=  0.0153  with  15-2=13  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  1.21144  1  1.2114  5162 

Residual  0.003051  13  0.000235 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constant  -1.81778  0.0145  -125 

logVn  1.79686  0.0250  71.8 

SPF6"=0.015213*(V/1000r  1.79686 
SPF56/100  =  0.361497*(V/1000)'^1. 79686 


SPF  7" 


SPF  7"  collected  18  Jan90  analyzed  again  2  April  1990 

SP7"L14/100=  0.215595*(V/1000)''1.84332 
SPL7"12/100  =0.224100*  (V/1000)''1.8473 
SPL7"23/100  =  0.214152*(V/1000)'^1.80784 

■:)<********  SP14  ***************** 

Summary  statistics  for  errPrSP14 

cases  selected  according  To  selV 
NumNumeric  =  12 
Mean  =  0.00009 
Standard  Deviation  =  0.01073 

Depencent  variable  is:  logcmSP14 

cases  selected  according  To  selV 
1 3  total  cases  of  which  1  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0058  with  12  -  2  =  10  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  1.92434  1  1.9243  58046 

Residual  0.000332  10  0.000033 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 
Constant  -1.04897  0.0041  -256 

logVn  1.84332  0.0077  241 
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SPF  7"  cont 

5PP  "7"  L12  ********* 

Summary  statistics  for  errPrSP14 

cases  selected  according  To  selV 
NumNumeric  =  12 
Mean  =  -0.00042 
Standard  Deviation  =  0.00657 


Dependent  variable  is:  logcmSP12 

cases  selected  according  To  selV 
1 3  total  cases  of  which  1  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0062  with  12  -  2  =  10  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  1.93265  1  1.9326  49595 

Residual  0.000390  10  0.000039 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 
Constant  -1.28100  0.0044  -289 

logVn  1.84730  0.0083  223 

SPF  7"  L12  =  0.05236  *  (V/ 1000)^  1.847 3 
SPF  7"  L12/100  =  0.05236*100/(280.375/12) 

SPF  7"  L12/100  =0.224100*  (V/1000)'' 1.8473 

'**4<********  SP7  23  no  couplers  ******** 

Summary  statistics  for  errPr23 

cases  selected  according  To  selV 
NumNumeric  =  12 
Mean  =  -0.00082 
Standard  Deviation  =  0.00566 

Depiendent  variable  is:  logcmSP23 

cases  selected  according  To  selV 
13  total  cases  of  which  1  are  missing 
r2  =  99.9%  R2(adjusted)  =  99.9% 
s  =  0.0132  with  12  -  2  =  10  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  1.85097  1  1.8510  10678 

Residual  0.001733  10  0.000173 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 
Constant  -1.74147  0.0094  -186 

logVn  1.80784  0.0175  103 

PrSP23  =  0.01 8 136*(V/1000)^1. 80784 

PrSP23perlOO  =  0.018136*100/(101.625/12)*(V/1000)^1.80784 
PrSP23perlOO  =  0.2 14 152*(V/1000)-' 1.80784 
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SPF  7"  cont 

gp  (;OUplCr  ********** 

Dependent  variable  is:  estSPcoupler 

cases  selected  according  To  selV>20()0 
13  total  cases  of  which  5  are  missing 
=  99.4%  R2(adjusted)  =  99.3% 
s  =  0.0009  with  8-2  =  6  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-rau  "> 

Regression  0.000909  1  0.000909  1039 

Residual  0.000005  6  0.000001 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constant  -0.004680  0.0007  -6.61 

uVP  0.019560  0.0006  32.2 

Dependent  variable  is;  estSPcoupler 

cases  selected  according  To  selV>20()0 
1 3  total  cases  of  which  5  are  missing 
r2  =  98.5%  R2(adjusted)  =  98.3% 
s  =  0.0025  with  8-1=7  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  0.002795  1  0.002795  450 

Residual  0.000043  7  0.000006 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

uVP  0.016015  0.0008  21.2 
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SPF  10”  analyzed  4Feb90 

note:  SP12  here  is  the  longest  run  upstream  of  the  station,  same  as  SP14  on  other  diameters 
case  46:  2  Jan90  -  appears  to  be  a  typo.  Replaced  with  transducer  value 
SP12/100  =  0.144993*(V/1000)''l. 88942 

'**********  L12  corresponds  to  SP13  for  Jim  and  L14  for  other  ducts  ** 

Summary  statistics  for  errFh-cSPI2 

NumNumeric  =  20 

Mean  =  0.00041 

Standard  Deviation  =  0.01 105 


Dependent  variable  is:  logSP12  doesn't  include  12/01/89 

35  total  cases  of  which  15  are  missing 
R2  =  99.8%  R2(adjusted)  =  99.8% 
s=  0.0159  with  20-2=18  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F- ratio 

Regression  1.97527  1  1.9753  7818 

Residual  0.004548  18  0.000253 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constant  -1.22126  0.0102  -120 

log(V/1000)  1.88942  0.0214  88.4 

SP12/100  =  10^-1.22126*100/(497.25/12)*(V/1000)^l. 88942 
SP12/100  =  0.144993*(V/1000)^1.88942 

about  86%  of  predicted  by  Colebrook  for  e=0.005" 


SPF  4" 


SPF4"  data  16Jan90 
analyzed  14  Feb  90 

1200<V<6000 

Summary  statistics  for  Prerr 

NumNumeric  =  13 
Mean  =  0.00211 
Standard  Deviation  =  0.02457 


Dependent  variable  is:  LogSP14 

cases  selected  according  To  selV 

13  total  cases  of  which  1  are  missing 

R2  =  100.0%  R2(adjusted)  =  100.0% 

s  =  0.0041  with  12  -  2  =  10  degrees  of  freedom 


Source  Sum  of  Squares 

Regression  1.54516 

Residual  0.000164 


df  Mean  Square 

1  1.5452 

10  0.000016 


Variable 

Constant 

logV 


Coefficient  s.e.  of  Coeff  t-ratio 
-0.732919  0.0032  -228 

1.85094  0.0060  307 


SPF4"  =  1 0^-0.7329 19*(V/1000)'' 1.85094 
SPF4"  =  0.1850  *  (V/1000)^1.851 


Dejjendent  variable  is:  cSP14 

cases  selected  according  To  1250<selV<5500 

13  total  cases  of  which  2  are  missing 

r2  =  100.0%  R2(adjusted)  =  100.0% 

s  =  0.0237  with  11-2  =  9  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square 

Regression  15.4811  1  15.48  27638 

Residual  0.005041  9  0.000560 


Variable 

Constant 

PrSP14 


Coefficient  s.e.  of  Coeff  t-ratio 

0.001316  0.0131  0.101 

1.00074  0.0060  166 


F-ratio 

93958 


F-ratio 


SPF  5" 


130 


analyzed  13  Feb90 
collected  12  Jan90 
SPFS"  =  .3268  (V/lOOOr  1.85843 

Summary  statistics  for  errPrl4 

NumNumeric  =  12 
Mean  = -0.00179 
Standard  Deviation  =  0.00791 

No  Ff  value  erred  by  more  than  1% 

Dependent  variable  is:  logcSP14 

75  total  cases  of  which  63  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0038  with  12  -  2  =  10  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  2.10272  1  2.1027  142687 

Residual  0.000147  10  0.000015 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constant  -0.854544  0.0022  -384 

LogV  1.85843  0.0049  378 

SPF5"  =  10^-0.854544  (V/1000)''l.85843  *  100/(513.25/12) 
SPF5"  =  .3268  (V/1000)''1.85843 


Dependent  variable  is:  PredSPF 

75  total  cases  of  which  63  are  missing 
R2  =  100.0%  R2(adjusted)  =  100.0% 
s  =  0.0058  with  12  -  2  =  10  degrees  of  freedom 

Source  Sum  of  Squares  df  Mean  Square  F-ratio 

Regression  8.38512  1  8.3851  246134 

Residual  0.000341  10  0.000034 

Variable  Coefficient  s.e.  of  Coeff  t-ratio 

Constant  0.005201  0.0028  1.88 

cSP14  0.993536  0.0020  496 


Appendix  N 
Selected  Raw  Data 


JUNCTION  DATA 

.  ! 

Setup  Va/Vs 

Vb/Vs 

Va/Vb 

LPmass 

LPJmass 

1 

0.8 

0.9 

0.9 

2220.8 

829.5 

2 

1.9 

0.6 

3.1 

1848.7 

1044.2 

3 

2.1 

0.6 

3.7 

1100.8 

651.8 

4 

2.3 

2.3 

1 

2063.6 

1 137.2 

5 

1 

2.2 

0.4 

1364.2 

769.3 

6 

0.4 

0.4 

1 

1854.1 

764.9 

7 

0.4 

0.4 

1 

442.9 

170.1 

8 

1 

1 

1 

1810.8 

717.3 

9 

0.3 

0.3 

1 

819.3 

335.1 

1  0 

0.6 

0.6 

1 

1411.9 

531.1 

1  1 

0.8 

0.8 

1.1 

1  442 

517.3 

1  2 

0.7 

0.7 

1 

1304.7 

468.8 

1  3 

0.5 

0.5 

1 

1277.6 

495.7 

1  4 

2.1 

2 

.1 

901.6 

466.6 

1  5 

1.7 

1.7 

1 

908.4 

430.7 

1  6 

1.9 

1.8 

1 

904.7 

441.2 

1  7 

1.5 

1.5 

1 

1083.8 

501.8 

1  8 

1.3 

1.3 

1 

1 120.1 

442 

1  9 

1.2 

1.2 

1 

1167.8 

449.9 

20 

1 

1.1 

0.9 

1224.2 

431 

21 

1 

1.2 

0.8 

1045.4 

386.4 

22 

1 

1.4 

0.7 

943.2 

375.5 

23 

1 

1.9 

0.5 

754 

362.9 

24 

0.9 

0.9 

1778.9 

659.4 

25 

0.8 

0.9 

0.8 

1714.6 

^2 

26 

0.7 

0.7 

1635.2 

561.4 

27 

0.6 

1 

0.6 

1529.4 

566.8 

28 

6^ 

1 

0.5 

1515.8 

531.8 

29 

0.9 

2.5 

0.4 

1065.3 

610.8 

30 

0.4 

1 

0.4 

1821 

676.7 

31 

2 

0.9 

2.1 

1109.6 

549.6 

32 

0.4 

1 

0.4 

1895.1 

727.8 

33 

1.5 

1.5 

1 

1055.2 

443.2 

34 

1 

0.9 

1 

1322.7 

464.2 

35 

0 

0 

680.81 

273.6 

32 


DUCT  A 

lSetuf:#VPpts  VPavg 

Q 

TP 

SPF 

Power 

LPF 

1 

30 

0.828 

500.1 

-3.445 

0.124 

-1  722.8 

61  .8 

2 

30 

1.878 

760 

-7.606 

0.264 

-5781 

200.9 

3 

30 

1.367 

647.9 

-4.833 

0.197 

-3131.1 

127.5 

4 

30 

1.381 

651.6 

-4.756 

0.199 

-3098.8 

129.5 

5 

30 

0.304 

305.9. 

-6.491 

0.049 

-1986 

14.9 

6 

30 

0.268 

287.6 

-6.321 

0.043 

-1818 

12.5 

7 

30 

0.106 

179.4 

-2.443 

0.018 

-438.2 

3.3 

8 

30 

0.978 

545; 

-3.317 

0.144 

-1807.5 

78.6 

9 

30 

0.13 

198.9; 

-5.187 

0.022 

-1031  .4 

4.4 

1  0 

30 

0.439 

365.7: 

-3.89 

0.069 

-1422.3 

25.1 

1  1 

30 

0.698 

459.91 

-2.579 

0.105 

-1186 

48.5 

1  2 

30 

0.557 

410.8 

-2.103 

0.085 

-863.9 

35.1 

1  3 

30 

0.302 

300. 8i 

-4.09 

0.049 

-1230.2 

14.6 

1  4 

30 

0.822 

496.1 

-2.963 

0.123 

-1469.8 

61 

1  5 

30 

0.795 

487.7; 

-2.896 

0.119 

-1412.5 

58.1 

1  6 

30 

0.81 

492.3: 

-2.921 

0.121 

-1437.7 

59.7 

1  7 

30 

0.847 

502. Il 

-3.124 

0.126 

-1568.4 

63.5 

1  8 

30 

0.827 

497.71 

-3.037 

0.124 

-151 1 .7 

61 .5 

1  9 

30 

0.821 

497.51 

-2.94 

0.123 

-1462.7 

61.1 

20 

30 

0.711 

4631 

-2.972 

0.107 

-1375.7' 

49.7 

21 

30 

0.554 

409.11 

-3.381 

0.085 

■1383.3 

34.8 

22 

30 

0.447' 

365. 5i 

-3.779 

0.07 

-1381 

25.5 

23 

30 

0.266 

282. 6i 

-4.305 

0.043 

-1216.5 

12.2 

24 

30 

0.746 

478.61 

-3.517 

0.112 

-1683.3 

53.7 

25 

30 

0.614 

434.41 

-3.752 

0.094 

-1630 

40.7 

26 

30 

0.496 

385.41 

-4.05 

0.077 

-1560.8 

29.7 

27 

30 

0.344 

325.51 

-4.253 

0.055 

-1384.3 

17.8 

28 

0.289 

293. 5i 

-4.576 

0.047 

-1343.1 

13.7 

29 

30 

0.195 

242.41 

-6.409 

0.032 

-1553.3 

7.8 

30 

30 

0.174 

228.41 

-5.629 

0.029 

-1285.6 

6.6 

31 

_ ^ 

1.385 

649.61 

-4.983 

0.199 

-3236.9 

129.5 

32 

30 

0.177 

232.5! 

-5.722 

0.029 

-1330.3 

6.8 

33 

0.84 

504. 9i 

-3.108 

0.125 

-1569.3 

63.3 

34 

30 

0.752 

477.51 

-2.845 

0.113 

-1358.5 

54 

35 

X 

0 

o' 

-3.701 

0 

0 

0 

33 


3^ 


DUCTB 

Setup#VPi 


its  VPavc 


0.989 


546.2  -3.151 


0.146 


Power 

-1721 


2 

30 

0.194! 

244.8 

-9.561! 

0.032 

-2340.1 

7.9 

3 

30 

0.101 

175.7 

-6.312! 

0.017 

-1  109.2! 

3.1 

4 

30 

1.44 

665.3' 

-4.615; 

0.207 

-3070.7 

137.5 

1.531 

0.272 

^111! 

0.8961 

0.141 

0.4291 

0.61  8T 
0.572! 

0.307 
0.76 
0.746| 
0.7481 
0.801 1 
0.782~ 
0.763  ~ 
0.812~ 
0.87" 
0.941  ' 

1.005' 

0.941  ' 
0.958' 


1.019 

0.995 

1.053 

1.453 


1.268 

0.302  ~ 

1.293~ 

0.819 

0.736 

0.742 


686.3: 

289.6' 
184: 
521. 8i 
207' 


-4.976 

-6.232 

-2.404 

-3.364 

-5.153; 


361.3!  -3.8781 


433.21 

416.9' 

303.3! 

476.9! 

472.4I 

472. 5I 

488.6! 

484.5! 

479.8  ~ 

495.1  ' 

512.4 

530.7' 

548.7! 

537.8' 

542.61 


66 


616.6 

303.6 

628.1 

498.4 

472.4 
474.31 


-2.628: 

-2.049! 

-4.04i 

-2.997' 

-2.925j 

-2.964j 

-3.159! 

-3.052! 

-2.979 

-2.81 

-2.95 

-3.131 

-3.357 

-3.212 

-3.277 


-4.168 
-6.286| 
-4.21 1! 
-3.104i 
-2.8281 
-2.6831 


0.219 

0.044 

0.019 

0.133 

0.024 

0.067 

0.094 

0.088' 

0.0491 

0.114i 

0.112! 

0.113j 
0.12 
0.1 17 
0.115 
0.121 
0.129 
0.139 
0.148 
0.139 
0.141 
0.15 
0.147 
0.155 
0.209 


0.184 

0.048 

0.187' 

0.122 

0.111 

0.112 


-3415.4 

-1804.8 

-442.3 
-1755.6: 
-1066.4 
-1401  .4 
-1  138.6 
-854 
-1225.3 
-1429.5 
-1381 .8' 
-1400.7' 
-1543.3, 
-1478.6! 
-1429.5 
-1391 .1 
-151 1 .7! 
-1661. 9j 
-1841 .61 
-1727.5! 

-17781 

-1852.9! 

-1875.8! 

-1986' 
-31  74.7' 


-2570.4! 
-1908. si 
-2644.9! 
-1546.9! 
-1335.81 
-1272.4! 


150.1 

12.7 
3.5 

69.5 
4.9 

24.2 

40.8 

36.6 

14.9 
54.5 

53.1 

53.2 

58.7 

56.9 
[ _ 55 

60.1 

66.3 
_ 1± 

81 .3 

74.9 

76.7 

82.8 
81.1 
86.8 

137.9 


113.4 
14.7 

117.5 

_ 6JL 

52.3 

52.9 


DUCT  S  ; 

Setup  #\/PptsVPavc 


1 .331  i 
0.528I 
0.321  ~ 
0.267~ 
0.32" 
1.984I 
0.734! 


1  242! 
789.6 
614.91 
561.3! 
615.6! 
1532.6! 
926.7! 


0.903!  1026.5; 


1.168! 

1.306! 

1.025! 

1.026! 


1.385 


0.1881 


1  169i 
1234.9; 
1091.8: 
1092.3! 


1261.2 


465.7 


-3.112 
-9.25! 
-6.126i 
-6.09! 
-6.5051 
-4.672 
-1 .835 
-3.549 
-4.155! 
-3.1051 
-2.393 
-1 .762! 


SPF 

0.074 

0.031 

0.02 

0.017 

0.02i 
0.106 
0.043; 
0.052 
_ 0^065i 

_ g. 

_ 0. 

0. 


Power 
-3865~^ 
-7304 
-3767 
-341  9| 
-4005 
-71 61 1 
-1701! 
-3643| 
-4857 


0.566 


0.7061 


0.574! 


809.3 


905.1 


816. 1i 


-3.3881 


-3.1571 


-3.499! 


1.048 


1.066 


1.131  [ 
1.096j 
1.16! 


0.231 

1.392' 

0.336' 

1.393' 

0.381  ' 

0.817' 

0.799' 


1112.7 


1122.2 


1140.3 
1137.61 
1 152.6 


516.5 

1267.4' 

628' 

1278.8" 

666.3' 

975.7' 

964.61 


-3.6061 


01 
-6.415 
-4.489~ 
-6.34' 
-4.58' 
•3.752' 
-2.9681 
-2.932! 


0.041!  -28 
0.034!  -28 


22i 

20 

0.4691 

734.41 

-3.8661 

0.028'  -: 

20 

0.2771 

564.7 

-4.357! 

MitMM 

0.059 

0.06" 

0.064' 

0.062! 


65! 


0.015 

0.077' 

0.021  ' 

0.077 

0.023 

0.047 

0.046! 


•37981 

•3900 

•4075i' 

•41021 


3313 
5690' 
3981  ' 
5856' 
2500' 
89 


•28281 


136 


DUCTM  ^  ^  ^ 

1  1 

1 

Setup  #VPpt^Q 

Qmass 

VPavgmasiTPmass 

SPFmassiPowermasdLPFmass 

1  i  301  2308.7 

2288.2 

1.0851  -4.165 

0.5061  -9530.3,  1  158.2 

2  30|  1852.4 

1794.4 

0.6541  -9.626 

0.318!  -17273.4'  570.9 

3  30|  1462.6 

1438.6 

0.421  -6.331 

0.213 

-9108.3 

306.3 

4 

30 

1871.3 

1878.2 

0.717 

-6.204 

0.346 

-1 1651 .7 

650 

5 

30 

1599.1 

1607.9 

0.524 

-6.698 

0.26 

-10770.5 

417.7 

6 

30 

2126.8 

2109.8 

0.903 

-5.99 

0.427 

-12637.7 

901.1 

7i  30 

1279.1 

1290 

0.341 

-2.344 

0.176 

-3024.2 

226.6 

CD 

CO 

o 

2061.2 

2093.3' 

0.9 

-4.308 

0.426 

-9017.4 

892.5 

9!  30 

1570.3 

1574.9 

0.509 

-4.937 

0.253 

-7774.4 

398.6 

1  Oj  30 

1984.2 

1961.9 

0.79 

-4.113 

0.378 

-8069.9 

742.1 

1  1 '  30 

2006.4 

1985 

0.814 

-3.214 

0.389 

-6379.6 

772.2 

1  2 

30 

1945.2 

1920 

0.759 

-2.577 

0.365 

-4947.7 

700.9 

1  3 

30 

1892.9 

1865.3 

0.728 

-4.089 

0.352 

-7626.5 

655.9 

1  4 

30 

1441.8 

1438.7 

0.432 

-3.839 

0.218 

-5523.7 

313.9 

1  5 

30 

1507.4 

1506.8 

0.473 

-3.741 

0.237 

-5636.5 

357.5 

1  6 

30 

1457.6 

1484.7 

0.46 

-3.783 

0.231 

-5616.4 

342.8 

1  7 

30 

1634.7 

1626i  0.554 

-4.022 

0.274 

-6539.7 

446 

1  8 

30 

1751.1 

1739.6 

0.631 

-3.91 

0.309 

-6802.5}  537.1 

1  9 

30 

1800.3 

1786.6 

0.661 

-3.807 

0.322 

-6801.51  574.7 

20 

30 

1879.2 

1863.2 

0.718 

-3.676 

0.347 

-6848.8 

646.2 

21 

30 

1743.6 

1737.5 

0.624 

-3.91 1 

0.305 

-6796 

530.1 

22 

30 

1635.3 

1630.5 

0.555 

-4.186 

0.274 

-6825.4 

447.5 

23 

30 

1391.5 

1396 

0.406 

-4.494 

0.206 

-6272.9 

287.8 

24 

30 

2152.2 

2129.1 

0.92 

-4.221 

0.435 

-8987.8 

925.2 

25 

30 

2094.1 

2099.2 

0.897 

-4.298 

0.425 

-9022.7 

891.7 

26 

30 

2081.2 

2077.6 

0.902 

-4.392 

0.428 

-9123.9 

888.9 

27 

30 

2012.1 

2016.2 

0.826 

-4.41 

0.394 

-8891 .6 

793.7 

28 

30 

2004 

2006.9 

0.844 

-4.597 

0.403 

-9225.5 

808.5 

29 

30 

1407.4 

1419.7 

0.419 

-6.414 

0.212 

-9106.6 

301.1 

30 

30 

2109.6 

2112.4 

0.928 

-5.381 

0.439 

-1  1366.7 

926.9 

31 

30 

1595 

1581.2 

0.512 

-6.473 

0.255 

-10235.9 

402.8 

32 

30 

2126.3 

2139.4 

0.936 

-5.481 

0.442 

-1  1726.5 

944.8 

33 

30 

1665.8 

1669.6 

0.574 

-3.996 

-6671 .4 

472.1 

34 

30 

1939.2 

1925.6 

0.763 

-3.59 

0.367 

-6912.6 

35 

30 

1431.7 

1438.9 

0.426 

-3.323 

0.215 

-4781 

309.8 

UNIVERSITY  OF  WASHINGTON 

SEATTLE,  WASHINGTON  98195 
13  June  1990 


Department  of  Environmental  Health  SC-34 

School  of  Public  Health  and  Community  Medicine 


Captain  William  Edmonson 
Program  Manager 
AFTT/CIMI 

Wright- Patterson  Air  Force  Base,  OH  45433-6588 
Dear  Captain  Edmonson: 

This  is  to  inform  you  that  James  P.  Curran  has  successfully 
completed  the  requirements  for  the  Master  of  Science  degree  in  the 
Department  of  Environmental  Health,  School  of  Public  Health  and 
Community  Medicine,  at  the  University  of  Washington.  He  will  be 
awarded  the  degree  in  June  1990. 

Si^rely, 

Jack  B^Hatlen 
Associ^  Professor  and 
Graduate  Program  Coordinator 


JBH:mlw 


